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ABSTRACT 

Using semi-analytic models of galaxy formation, we investigate galaxy properties such 
as the Tully-Fisher relation, the B and K-band luminosity functions, cold gas con- 
tents, sizes, metallicities, and colours, and compare our results with observations of 
local galaxies. We investigate several different recipes for star formation and supernova 
feedback, including choices that are similar to the treatment in Kauffmann, White & 
Guiderdoni (1993) and Cole et al. (1994) as well as some new recipes. We obtain good 
agreement with all of the key local observations mentioned above. In particular, in 
our best models, we simultaneously produce good agreement with both the observed 
B and K-band luminosity functions and the I-band Tully-Fisher relation. Improved 
cooling and supernova feedback modelling, inclusion of dust extinction, and an im- 
proved Press-Schechter model all contribute to this success. We present results for 
several variants of the CDM family of cosmologies, and find that models with values 
of £1 — 0.3-0.5 give the best agreement with observations. 
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1 INTRODUCTION 

Over the past decade and a half, a great deal of progress to- 
wards a qualitative understanding of galaxy properties has 
been made within the framework of the Cold Dark Mat- 
ter (CDM) picture of structure formation (e.g., Blumentha 
|et al. 1984 ). However, N-body simulations with gas hydro 



dynamics still have difficulty r eproducing the o bserved prop- 
erties of galaxies in detail (cf. [gteinmetz 1997 ). It is appar- 
ent that there must be additional physics that needs to be 
included in order to obtain realistic galaxies in the CDM 
framework. It is likely that many processes (e.g. cooling, 
star formation, supernova feedback, etc.) form a complicated 
feedback loop. It is not computationally feasible to include 
realistiq physics over the required dynamic range in N-body 
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by the unknowns associated with star formation, supernova 
feedback, the stellar initial mass function, metallicity yield, 
dust extinction, etc. However it is not only a question of 
computational efficiency: the macroscopic picture afforded 
by the semi-analytic method provides an important level of 
understanding that would be difficult to achieve by running 
an N-body simulation, even if we had an arbitrarily large 
and fast computer. 

The semi -analytic approach to galaxy formation was 
formulated in White & Frcnk (1991), but this approach was 
not Monte-Carlo based and thus could only predict average 
quantities. The Monte-Carlo approach was primarily devel- 
oped independently by two main groups, which we shall refer 
to as the "Munich" group (Kauffmann, White, & Guider- 



doni 



1993j 



mann 1995 



Kauffmannj Guiderdoni 



Kauffmann 1996a; 



White 199 



Kauffmann 1996 



mann, Nusser, & Steinmetz 1997; Kauffmann fc Chariot 



Se mi-analytic models (SAMs) of galaxy formation are 1996) and t he "Durham" group ( Cole et al. 1994; Hey] 



embedc ed within the framework of a CDM-like initial power et al. 1995]; Baugh, Cole, & Frenk 1996a 



Kauff- 



Kauff- 



spectru m and the theory of the growth and collapse of fluc- 



Frenk 19961; Baugh et al. 1997), because the major- 



Baugh, Cole 



tuations through gravitational instability. They include a 
simplified yet physical treatment of gas cooling, star forma- 
tion, supernova feedback, and galaxy merging. The Monte- 
Carlo approach enables us to study individual objects or 
global quantities. Many realizations can be run in a mod- 
erate amount of time on a workstation. Thus SAMs are an 
efficient way of exploring the large parameter space occupied 



ity of the members of these groups are associated with the 
Max-Planck-Institut fur Astrophysik in Garching, near Mu- 
nich, Germany, and the University of Durham, U.K., re- 
spectively. Similar mode ls have also been investigated by 
Lacey fc Silk (199l| ) and |Lacey et al. (1993| ). This work has 
shown that it is possible to reproduce, at least qualitatively, 
many fundamental observations in the simple framework of 
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SAMs. These include the galaxy luminosity function (LF), 
the Tully- Fisher relation (TFR) , the morphology- density re- 
lation, cold gas content as a function of luminosity and en- 
vironment, and trends of galaxy colour with morphology 
and environment. However, some unsolved puzzles remain. A 
fundamental discrepancy has been the inability of the mod- 
els to simultaneously reproduce the observed Tully-Fisher 
relation and the B-band luminosity function in any CDM- 



type cosmology ( Kauffmann, W hite, & Guiderdoni 1993 



Cole et al. 1994 |Heyl et al. 1995| ). Another problem, thought 



to be generic to the hierarchical structure formation sce- 
nario, is the tendency of larger (more luminous) galaxies to 
have bluer colours than smaller (less luminous) ones, in con- 
trast to the observed trend. We shall discuss these and other 
problems in detail in this papeiQ. 

This paper has several goals. We describe the ingredi- 
ents of our models and show that they reproduce fundamen- 
tal observations of the local universe. We repeat the calcu- 
lations of several quantities that have been studied before 
using SAMs, and one might wonder why this is worthwhile. 
First, this will serve as a reference point for future papers in 
which we will use these models to study new problems. Sec- 
ond, the previous studies have been spread out over several 
years with different quantities being presented in different 
papers. Over this time the models themselves have evolved. 
We therefore think it will be useful to have all of these results 
presented in the same place in a homogeneous manner. In 
addition, the two main groups have not always studied the 
same quantities, and when they have, they have not always 
presented their results in a way that is directly comparable. 
This makes it difficult for the non-expert to judge just how 
different these two approaches really are. Moreover, because 
the models differ in so many details, it is impossible to deter- 
mine which particular ingredients are responsible for certain 
differences in the results. Two of the important differences 
that we are particularly interested in are the parameteri- 
zation of star formation and supernova feedback. We shall 
investigate the results of varying these recipes while keeping 
the other ingredients fixed. We also include some physical 
effects that have previously been neglected, and show that 
some of the problems that have plagued previous models can 
be alleviated. We investigate the importance of the under- 
lying cosmology by examining the same quantities in a wide 
range of different cosmologies, spanning currently popular 
variants of the Cold Dark Matter (CDM) family of models. 

The structure of the paper is as follows. In Section ^ 
we describe the basic physical ingredients of the models 
and briefly summarize the SAM approach. In Section |^, we 
summarize the model parameters and describe how we set 
them. In Section we illustrate the effects of varying the 
free parameters and the star formation and supernova feed- 
back recipes, using the properties of galaxies within a "Local 
Group" sized halo as an illustration. In Section [|, we present 
the results of our models for fundamental global quanti- 
ties and galaxy properties, illustrate the effects of different 
choices of star formation and supernova feedback recipes 



Naturally these groups have continued to modify and improve 
their models. In this paper, when we make general statements 
about the published Munich and Durham models, we refer to 
work that was published before February 1998 



Table 1. Parameters of Cosmological Models. From left to right, 
the tabulated quantities are: the matter density, the density in 
the form of a cosmological constant in units of the critical density, 
the Hubble parameter, the baryon density in units of the critical 
density, the age of the universe in Gyr, the slope of the primordial 
power spectrum, and the linear rms mass variance on a scale of 
Sh- 1 Mpc. 



Model 


n 


H A 


h 


U b 


to 


« 


0"8 


SCDM 


1.0 


0.0 


0.5 


0.072 


13.0 


1.0 


0.67 


tCDM 


1.0 


0.0 


0.5 


0.072 


13.0 


1.0 


0.60 


ACDM.5 


0.5 


0.5 


0.6 


0.050 


13.5 


0.9 


0.87 


OCDM.5 


0.5 


0.0 


0.6 


0.05 


12.3 


1.0 


0.85 


ACDM.3 


0.3 


0.7 


0.7 


0.037 


13.5 


1.0 


1.0 


OCDM.3 


0.3 


0.0 


0.7 


0.037 


11.3 


1.0 


0.85 



on these quantities, and compare our results with previous 
work. We summarize and discuss our results in Section ^. 



2 BASIC INGREDIENTS 

In this section we summarize the simplified but physical 
treatments of the basic physics used in our SAMs. This 
includes the growth of structure in the dark matter com- 
ponent, shock heating and radiative cooling of hot gas in 
virialized dark matter halos, the formation of stars from 
the cooled gas, the reheating of cold gas by supernova feed- 
back, the evolution of the stellar populations, and mergers 
of galaxies within the dark matter halos. There are many 
assumptions implicit in this modelling and in addition to de- 
scribing the choices we have adopted in our fiducial models, 
we also remark upon some relevant details of the assump- 
tions made in previously published work. 

Our models have been developed independently, but 



very much in the spirit of Kauffmann, White, fc Guider^ 
doni (1994 hereafter KWG9377 |Cole et al. (1994 hereafter 



CAFNZ94), and subsequent work by these groups. We refer 
the reader to this literature for a more detailed introduc- 
tion to the SAM approach, which here is summarized rather 
briefly. A more detailed review of the literature a nd descrip- 
tion o f an earlier version of our models is given in Somerville 



(1997) 



2.1 Cosmology 

Most of the previous SAM work has been in the context 
of standard cold dark matter (SCDM), Qq — 1, Ho = 
50 km s _1 Mpc -1 , as = 0.67. However, this model has now 
been discredited many times in many different ways. Partic- 
ularly relevant to our work is the problem that this model 
overproduces objects on galaxy scales relative to cluster 
scales. In addition the normalization ag = 0.67 is highly 
i nconsistent with th e COBE data, which requires as ~ 1.2 
(|G6rski et al. 1996) for this model. Many alternative vari- 



ants of CDM have been suggested. We have chosen illus- 
trative examples of popular variants of the CDM family 
of models, spanning the observationally plausible range of 
parameter space. We have retained the standard SCDM 
model for comparison with previous work, and consider one 
other model with fio = 1, the rCDM model of Efstathiou 



© 0000 RAS, MNRAS 000, 000-000 



Semi- Analytic Modelling of Galaxy Formation 3 



SCDM 

tCDM 

ACDM.5 

0CDM.5 

ACDM.3 



(>c'IAI.:s 




10 13 
(h-> M G 



Figure 1. The mass function From the imnroved Press-Schechter 



model p roposed by ( Shcth fc Tormcn 1999|) divided by the stan- 



dard Prcss-Schechter mass function. This "correction factor" is a 
function of redshift as well as halo mass, and here is shown for 

2 = 0. 



Bond, & White (1992). For our purposes, the properties of 
this model are very similar to other popular Q = 1 models 



LllLed CDM (n ~ 0.8) an d nrogefe wttfa an ad mix- 



hot daik .11, .in i (CIIDM, Pmaaek et al. 19 95 ). We 



such as 
ture of 

also consider open (Qa ~ 0) and flat Slo + = 1) models 
with Slo =■ 0.3 and Slo = 0.5. We have assumed a Hubble 
parameter h = 0.7 for the Qq = 0.3 models, h — 0.6 for 
the fio = 0.5 models, and h = 0.5 for the flo = 1 models 
(H = lOOftkm s _1 Mpc -1 ). For the tt = 0.5 fiat model, 
we have included a mild tilt (n = 0.9) to better simulta- 
neously fit the power on COBE and cluster scales. For all 
the models, in computing the power spectrum we have as- 
sumcd the baryon fraction implied by the observations of 
Tytler et al. (1999[), Q. h = 0.019 h -2 . All models assume 



T/S — (no contri bution from tensor mo des). We use the 
fitting functions of Bunn & White (1997), modified to ac- 
count for the presence of baryons using the prescription of 



Hu fc qugiyama (199t ), to obtain the linear power spectra 
and COBE iiuiiiiahzafiuiis. The iiumiahzafiuii u$ is luughly 
consistent with the z = cluster abundance and the COBE 
measurement except in the case of SCDM and the Qo = 0.3 
open m odel, for which we have uacd the cluatcr normaliza — 



tion. T he parameters of the cosmologies arc presented in 



Table 0. 



2.2 Dark Matter Merger Trees 



The extend e d Press-Schechter for malism (Bower 1991 



et al. 1991 



Bond 



Lacey & Cole 1993) provides us with an ex- 



pression for the probability that a halo of a given mass 
mo at redshift zo has a progenitor of mass mi at some 
larger redshift Zi. Several methods of creating Monte-Carlo 
realizations of the merging histories of dark matter halos 
( "merger trees") using this f ormalism h a ve been developed 



(Kauffmann fc White 1993; Cole 1991 jomerville fc Ko- 



latt 199£ ). Although the agreement of the Press-Schechter 
model with N-body simulations is in some ways surpris- 
ingly good given the simplifications involved, recent work 
has emphasized that there are non-negligible discrepancies. 



Several autho rs (|Gross et al. 1998|; |Somerville et al. 199£ 



Tormen 1998 



Ibzzi & Governato 1997) have now reported 



the same results using different N-body codes and different 
methods of identifying halos, indicating that the problems 



are un likely to be explained by numerical effects. Gross et al 
(1998[ ) showed that for a wide variety of CDM-type models 



the z — Press-Schechter mass function agrees well with 
simulations on mass scales M > 10 14 Mq, but on smaller 
scales the Press-Schechter theory over-predicts the number 
of halos by about a factor of 1.5 to 2. The precise factor 
varies somewhat depending on the cosmological model and 
power spectrum and the way the Press-Schechter model is 
implemented. However, this problem cannot be solved by 
adjusting the critical density for collapse, 8 c ,o- In addition, 
the Press-Schechter model predicts stronger evolution with 
redshift in the halo mass function than is observed in the 



simula tions (dross 1997; Somerville et al. 1998| ). [Tormen 
(1998| ) finds a similar behaviour when comparing the predic- 



tions of the extended Press-Schechter theory with the con- 
ditional mass function of cluster-sized halos in simulations. 



Somerville et al. (1998) also compared the extended Press- 
Schechter model with the results of dissipationless N-body 
simulations, and investigated how well the distribution of 
progenitor number and mass in the simulation s agrees with 
that produced by the merger-tree method of Somerville & 



Kolatt (1998). They found that the distributions of pro- 



genitor number and mass obtained in the merging trees, 
which have been deliberately constructed to reproduce the 
Press-Schechter model, are skewed towards larger numbers 
of smaller mass progenitors than are found in the simula- 
tions. This problem is endemic to any method based on the 
extended Press-Schechter model. However, the relative prop- 
erties of progenitors within a halo of a given mass are very 
similar in the merger trees and the simulations. This sug- 
gests that the merger trees should provide a fairly reliable 
framework for modelling galaxy formation, if the overall er- 
ror in the Press-Schechter mass function is corrected for. 
An improved version of the Press-Schechter model, which 
gives good agreement with simulations fo r a variety of cos- 



molog ies, has recently been proposed by Sheth & Tormcn 



(19991). The "correction factor", i.e. the mass function from 



the Sheth- Tormen model divided by the standard Press- 
Schechter mass function, is shown in Fig, jlj. 



In the merger-tree method of somerville & Kolatt 

( 199S| ) , used here, the merging history of a dark matter halo 



is constructed by sampling the paths of indivi dual particle 



traje c tories using the ex cursion set formalism (Bond et al 



1991; Lacey & Cole 1995). It does not require the imposition 



of a grid in mass or redshift, nor are merger events required 
to be binary. The redshifts of branching events (i.e. halo 
mergers) and the masses of the progenitor halos at each 
stage are chosen randomly using Monte-Carlo techniques, 
such that the overall distribution satisfies the average pre- 
dicted by the extended Press-Schechter theory. Thus when 
we subsequently refer to a "realization" we mean a particu- 
lar Monte-Carlo realization of the halo merging history. This 
is the most important stochastic ingredient in the models. 
In order to make the tree finite, it is necessary to impose 
a minimum mass m m i n . Although the contribution of mass 
from halos smaller than m m i n is included, we do not trace 
the merging history of halos with masses less than m m i n , 
but rather assume that this mass is accreted as a diffuse 
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Figure 2, Cooling radius " f as a junction of circular kb!oc= 

ity. The straight diagonal line shows the virial radius, which the 



cooling radius may not exceed. The curved lines show the cool- 
ing radius predicted by the literal static halo cooling model (see 
text), assuming that the hot gas has primordial, 0.3 solar, or solar 
metallicity. Open circles show the application of the static halo 
model within the merger trees, and crosses show the dynamic halo 
model (see text), assuming a fixed metallicity of 0.3 solar. Earlier 
conversion of gas from the hot to cold phase and reheating of hot 
gas by halo mergers results in a lower cooling efficiency for large 
halos in the dynamic halo model. 



component. Here we take m m i n to be equal to the mass cor- 
responding to a halo with a circular velocity of 40 km/s at 
the relevant redshift. We argue that galaxies are unlikely to 
form in halos smaller than thi s because the gas will be pho- 
toio nized and unable to c ool (Weinberg, Hernquist, & Katz 
19971; iFlarcado-Miro 19971). 



For the prediction of global quantities, we run a grid of 
halo masses (typically ~ 50 halos from 10 m m i n to V c — 1500 
km/s), and weight the results with the overall number den- 
sity for the appropriate mass and redshift, using the im- 
proved Press- Schecter model of 3heth & Tormen (1999). We 
run many such grids and average the results. 

2.3 Gas Cooling 

2.3.1 Cooling in Static Halos 

Gas cooling is modelled using an app roach similar to the one 
introduced by White & Frenk (1991). We assume that each 
newly formed halo, at the top level of the tree, contains 
pristine hot gas that has been shock heated to the virial 
temperature of the halo (T) and that the gas traces the 
dark matter. The rate of specific energy loss due to radiative 
cooling is given by the cooling function A(T). We can then 
derive an expression for the critical density which will enable 
the gas to cool within a timescale r coo i: 

3 [im p UbT 



Pcool 



2 xl r C ooiA(T) 



(1) 



where pm p is the mean molecular weight of the gas and \ e = 
n e /ntat is the number of electrons per particle. Assuming 
that the gas is fully ionized and has a helium fraction by 
mass Y = 0.25, 

k B T 



pcool 



= 3.52 x 10 7 - 



ooiA 23 (T)' 



(2) 



where ksT is in degrees Kelvin, r coo i is in Gyr, and A23(T) = 
A(T)/(10 -23 ergs s _1 cm 3 ). The virial temperature is ap- 
proximated as ksT — 71.8<7 2 ; r , where is the virial ve- 
locity dispersion of the halo. If we assume a form for the 
gas density profile p g (r), we can now invert this expression 
to obtain the "cooling radius" , defined as the radius within 
which the gas has had time to cool within the timescale T CO oi . 
For the simplified choice of the singular isothermal sphere, 
this gives 



?~cool — 



P0 
pcool 



1/2 



(3) 



where po = /hot V - c 2 /(47rG), /hot is the hot gas fraction in the 
cooling front and V c = \/2o"vir is the circular velocity of the 
halo. 



We use the cooling function of Sutherland & Dopita 



(1993| ). The cooling function A(T) is also metallicity depen- 



dent. In this paper, we assume that the hot gas has an aver- 
age [Fe/H] = 0.3 Z@ at all redshifts and for all halo masses. 
This value is typical of the hot gas in clust ers from z — 
to z ~ 0.3 (Mushotzky & Loewenstein 1997). We will treat 



chemical enrichment in more detail and investigate the ef- 
fects on cooling in a future paper. In practice, we find that 
the results at z = are very similar regardless of whether 
we use the self-consistently modelled hot gas metallicity in 
the cooling function, or a fixed metallicity [Fe/H] = 0.3 Zq. 

We divide the time interval between halo mergers 
(branchings) into small time-steps. For a time-step At, 
the cooling radius increases by an amount Ar and we as- 
sume that the mass of gas that cools in this time-step is 
dm C ooi = 47rr- 2 ool p s (r C ooi) Ar. The cooling radius is not al- 
lowed to exceed the virial radius, and the amount of gas 
that can cool in a given timestep is not allowed to exceed 
the hot gas contained within the virial radius of the halo. 
For small halos, and at high redshift, the cooling is therefore 
effectively limited by the accretion rate. New hot gas is con- 
stantly accreted as the halo grows. When we construct the 
merging tree, we keep track of the amount of diffuse mass 
(i.e. halos below the minimal mass m m i n ) accreted at ev- 
ery branching, m acc . The mass of hot gas accreted between 
branchings is then / bar m acc , where / bar = fib/^o is the uni- 
versal baryon fraction. We assume that the mass accretion 
rate is constant over the time interval between branchings, 
which is what one would expect from the spherical collapse 
model (see Appendix). We also require that even if the gas 
is able to cool, it falls onto the disk at a rate given by the 
sound speed of the gas, c s = (5kBT/3pm p ) 1 ^ 2 ~ 1.3cr„, 
where a v is the 1-D velocity dispersion of the halo. Note 
that c 3 is approximately equal to the dynamical velocity of 
the halo, and that N-body simulations with hydrodynamics 
and cooling show that the radial infall velocity of cooling 
gas within the virial radius is gen erally close to this value 
(Evrard, Metzler, fc Navarro 199^). 
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2.3.2 Cooling and Heating in Merging Halos 

In the simplest version of this approach to modelling gas 
cooling in dark matter halos, we imagine that a halo of a 
given circular velocity V c , with a corresponding virial tem- 
perature T, forms at time t = and grows isothermally, 
gradually cooling at the rate given by the cooling function 
A(T) as described above. In this model, the cooling time 
Tcooi is the age of the Universe at the current redshift, and 
the gas fraction in the cooling front /hot is always equal to 
the universal baryon fraction /bar- We refer to this picture 
as the "static halo" cooling model, as it does not account 
for the dynamical effects of halo mergers. 

However, in the hierarchical framework of the merger 
trees, most halos are built up from merging halos that have 
experienced cooling, star formation, and feedback in previ- 
ous time-steps. This modifies the gas fraction /hot in the 
cooling front. The virial velocity and temperature change 
discontinuously following a merger, and merger events may 
shock heat the cooling gas. We have developed a "dynamic 
halo" cooling model that incorporates these effects in the 
following way. 

For top-level halos (halos with all progenitors smaller 
than the minimal mass m m in), the gas fraction /hot is as- 
sumed to be equal to the universal baryon fraction /bar , and 
the cooling time t coo i is the time elapsed since the initial col- 
lapse of the halo. Subsequently, when a halo forms from the 
merging of two or more halos larger than m m i n , we determine 
whether the mass of the largest progenitor mi comprises 
more than a fraction / re hoat of the post-merger mass mo. If 
so, the cooling radius and cooling time of the new halo are 
set equal to those of the largest progenitor. The gas fraction 
in the cooling front is taken to be /hot = mhot/nitot(> r coo i), 
where mhot is the sum of the hot gas masses of all the pro- 
genitors, and mtot(f > r coo i) is the total mass contained 
between the cooling radius and the virial radius of the halo 
(assuming an isothermal profile). If mi /mo < /reheat, we 
assume that the hot gas within all the progenitor halos is 
reheated to the virial temperature of the new halo, and the 
cooling radius and cooling time are reset to zero. The gas 
fraction in the cooling front is then /hot = mhot /mo- Note 
that /hot may in principle be larger than /b ar due to re- 
heating by supernova, but in general /hot < ./bar because of 
previous gas cooling and consumption. 

We can apply the main simplifying assumptions of the 
static halo cooling model within the merging trees, i.e. we 
always assume /hot = /bar and t coo i equal to the age of the 
Universe at any given time, and do not reheat the gas after 
any halo mergers. The results differ somewhat from the lit- 
eral static halo model because we do not allow the cooling 
rate to exceed the available supply of hot gas, or to exceed 
the sound speed constraint, and because the progenitor ha- 
los cool at different temperatures. Fig. ^| shows the cooling 
radius in the literal static cooling model, and in the appli- 
cation of the static cooling model within the merging trees. 
For low V c halos, the cooling is limited by the available col- 
lapsed gas supply (i.e. r coo i > r v i T ). For larger halos the 
results are similar to the prediction of the literal static halo 
model. However, the dynamic halo model (crosses) predicts 
significantly less cooling in large halos, due to the lower val- 
ues of /hot and the reheating by halo mergers. Note that 
the cooling model used by the Munich group more closely 



resembles the "static halo" model, and the cooling model 
used by the Durham group is more similar to our "dynamic 
halo" cooling model^ In this paper, we will show results for 
both cooling models. 



2.4 Disk Sizes 

To obtain a very rough estimate of the sizes of disks 
th at form in our models, w e adopt the general picture 
of Fall & Efstathiou (198C), in which the gas collapse 



is halted by angular momentum conservation. We define 
A_h to be the dimensionless spin parameter of the halo, 
Ah = J|-E| 1/2 G -1 M" 5/2 , where J is the angular momen- 
tum, M is the mass and E is the energy of the dark matter 
halo. We assume that the gas has the same specific angu- 
lar momentum as the dark matter, and collapses to form 
a disk with an exponential profile. For a singular isother- 
mal halo, the scale radius of the disk that forms is then 
r s — 1/ '\f2\uri, where rt is the radius before collapse (in 
our models, n = min(r coo i, r vir ). We neglect the modifica- 
tion of the inner profile of the dark matter due to the i nfall 



of the baryons, whic h will tend to lead to smaller disks ( Blu- 
Flores et al 



menthal et al. 1986 



199S 



1993 



Mo, Mao, & White 



The distribution of Ah found in N-body simulations 



3arnes & Efstathiou 1987; Kravtsov et al. 1998 



Lemson 



fa Kauffman 1997) is a rather broad log- normal with mean 



(Ah) = 0.05. It is likely that in order to obtain a realistic dis- 
tribution of galaxy sizes and surface brightnesses, we should 
consider a range of values of Ah as seen in the above simu- 
lations (tDalcanton, Spergel, fc Summers 1997 ; Mo, Mao, & 
White 1998[ ). However, it is not known how Ah is affected 
by mergers, so we do not know how to propagate this quan- 
tity through the merging trees. We should also use a more 
realistic halo profile than the singular isothermal sphere. We 
intend to address the modelling of disk sizes in more detail 
in the future. For the present, we use Ah = 0.05 for all halos. 



2.5 Star Formation 

The star formation recipes that we will consider in this paper 
are of the general form 



m co id 



(4) 



where m co id is the total mass of cold gas in the disk and we 
hide all of our ignorance in the efficiency factor r» . The sim- 
plest possible choice is to assume that r* = r° is constant. 
This would imply that once it is cold, gas is converted to 
stars with the same efficiency in disks of all sizes and at all 
redshifts. We shall refer to this recipe as SFR-C. 

Another choice is a power law, in which the star for- 
mation efficiency is a function of the circular velocity of the 
galaxy: 



"Icold 



(5) 



t No te that in earlier versions of our models (e.g. Somervillc 



1997), as in the Munich models, we prevented gas from cooling 



altogether in large halos by applying an arbitrary cutoff. We no 
longer apply this cutoff. 
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Figure 3. The star formation rate per unit mass of cold gas 
(star formation efficiency) as a function of circular velocity. The 
solid curved line shows the Durham star formation law (SFR-D), 
which has no explicit dependence on redshift, and the horizontal 
dashed lines show the Munich star formation law (eqn. SFR-M), 
for z = 5, 3, 1, from highest to lowest, respectively. 



T.(V C ) 



fV c \ 



(6) 



where m co id is the mass of cold gas in the disk, r* and a* 
are free parameters, and Vo = 300 km s _1 is an arbitrary 
normalization factor. This is equivalent to the approach used 
by the Durham group, and we will refer to this recipe as 
SFR-D. 

The other approach, used by the Munich group, as- 
sumes that the timescale for star formation is proportional 
to the dynamical time of the disk 



m co id 

tS Tdyn 



(7) 



Here r" is a dimensionless free parameter, and Tdyn is the 
dynamical time of the galactic disk, Tdyn = ^disk/14. Follow- 
ing KWG93, we take Tdisk to be equal to one tenth the virial 
radius of the dark matter halo, and V c to be the circular ve- 
locity of the halo at the virial radius. For satellite galaxies, 
the dynamical time remains fixed at the value it had when 
the galaxy was last a central galaxy. We will refer to this 
star formation recipe as SFR-M. 

It is worth noting the differences in these assumptions 
and the implications for the models. The dynamical time 
i"dyn at a given redshift is nearly independent of the galaxy 
circular velocity. This is because the spherical collapse model 
predicts that the virial radius scales like r v j r oc V c (see Ap- 
pendix). However, the virial radius of a halo with a given cir- 
cular velocity increases with time (r v i r oc (1 + z)~ 3 ^ 2 for an 
Einstein-de Sitter universe). This means that SFR-M is ap- 
proximately constant over circular velocity but has a higher 
efficiency at earlier times (higher redshift). In contrast, SFR- 
D has no explicit dependence on redshift but does depend 
fairly strongly on the galaxy circular velocity (a, = —1.5 in 



the fiducial Durham models), so that star formation is less 
efficient in halos with small V c . This is illustrated in Fig. ^. 
Because a "typical" halo at high redshift is less massive and 
hence has a smaller circular velocity in hierarchical models, 
this has the effect of delaying star formation until a later 
redshift, when larger disks start to form. SFR-M therefore 
leads to more early star formation. We will discuss this fur- 
ther in Section ^. 



2.6 Supernova Feedback 

2.6.1 Previous Feedback Recipes 

In the Munich and Durham models, the rate of reheating of 
cold gas is given by 



m rh = eg N (14/V )- a *m* 



(8) 



where eg N and a r h are free parameters, m* is the star for- 
mation rate, and Vo is a scaling factor chosen so that eg N 
is of order unity (Vo = 400 km s^ 1 for the Munich mod- 
els and 140 km s _1 for the Durham models). The Munich 
group assumes « rn = 2, whereas the fiducial models of the 
Durham group assume a considerably stronger dependence 
on circular velocity, a rn = 5.5. 

The "reheated" gas is removed from the cold gas reser- 
voir. An important issue is whether the reheated gas remains 
in the halo in the form of hot gas, where it will generally cool 
again on a short timescale, or is expelled from the potential 
well of the halo entirely. In the Munich models (previous to 
1998), all the reheated gas is retained in the halo (G. Kauff- 
mann, private communication). In the Durham models, all 
the reheated gas is ejected from the halo. This gas is then 
returned to the hot gas reservoir of the halo after the mass 
of the halo has doubled (Durham group, private commu- 
nication). We find that the results of the models are quite 
sensitive to whether the gas is retained in or ejected from 
the halo. We would therefore like to find a simple but phys- 
ical way of modelling the ejection of the reheated gas from 
the disk and the halo without introducing an additional free 
parameter. To this end, we have introduced the following 
modified treatment of supernova feedback. 



2.6.2 The Disk-Halo Feedback Model 

We assume that the mass profile of the disk is expo- 
nential. The potential energy of an exponential disk with 
scale radius r s and central surface density Sp is app roxi- 
mately W ~ -11.6GEgrf (|Binney fc Tremaine 198?|). We 



can then calculat e the rms escape velocity for the disk, 
(i>e SC , disk) 1//2 = v/ 4W/irid, where m d is the mass of the 
disk. Similarly, f or the halo, the rms escape velocity is 
(«esc, halo) 1/2 = \/-4W7m ha!o = VZVc, using the virial the- 
orem. As before, we have the free parameter eg N , which we 
interpret loosely as the fraction of the supernova energy 
transferred to the gas in the form of kinetic energy. The 
rate at which kinetic energy is transferred to the gas is now 
£sn = £sn -^sn ??sn rn* , where Esn = 10 51 ergs is the total 
(kinetic and thermal) energy per supernova, 7?sn is the num- 
ber of supernova per solar mass of stars (rjsN = 3.2 x 10~ 3 
for the Scalo IMF used here ( |Bruzual fc Chariot 199^ )), and 
m* is the star formation rate. Following the general argu- 
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ments of Dekel & Silk (1986), we now calculate the rate at 



which gas can escape from the disk: 

Tlrh.disk = ^/rh.disk-/— 5 T, 



(9) 



and from the halo (same expression with (v^ BC halo))- The 
factor / r h,disk is a fudge factor that we leave fixed to 
/rh.disk = 2. In general, dm r h,disk is then larger than 
dm r h,halo' The gas that can escape from the disk but not 
the halo is added to tlip hot gas in the halo Gas that ran 



escape Ifrom both the disk and the halo is removed from 
the halo entirely. Because of the uncertainties involved, we 
do not attempt to model the recollapse of this gas at later 
times, so this gas will never be re-incorporated into any halo 
and is "lost" forever. 



2.7 Chemical Evolution 

We trace chemical evolution by assuming a constant "effec- 
tive yield", or mean mass of metals produced per mass of 
stars. The value of the effective yield, y, is treated as a free 
parameter. We assume that newly produced metals are de- 
posited in the cold gas. Subsequently, the metals may be 
ejected from the disk and mixed with the hot halo gas, or 
ejected from the halo, in the same proportion as the reheated 
gas, according to the feedback model described above. The 
metallicity of each batch of new stars equals the metallicity 
of the cold gas at that moment. Note that because enriched 
gas may be ejected from the halo, and primordial gas is 
constantly being accreted by the halo, this approach is not 
equivalent to a standard "closed box" model of chemical 
evolution. Also note that although we track the metallic- 
ity of the hot gas by this procedure, in this paper we do 
not use this metallicity to compute the gas cooling rate (see 
Section 2.3). 



2.8 Galaxy Merging 

2.8.1 Dynamical Friction and Tidal Stripping 

When halos merge, we assume that the galaxies within them 
remain distinct for some time. In this way we eventually end 
up with many galaxies within a common dark matter halo, 
as in groups and clusters. The central galaxy of the largest 
progenitor halo becomes the new central galaxy and the 
other galaxies become "satellites" . Following a halo merger 
event, we assume that the satellites of the largest progeni- 
tor halo remain undisturbed and place the central galaxies 
of the other progenitors at a distance / mr g?" v i r from the cen- 
tral galaxy, where / mr g is a free parameter and r v i r is the 
virial radius of the new parent halo. Satellites of the other 
progenitors are distributed randomly around their previous 
central galaxy, preserving their relative distance from that 
galaxy. All the satellites lose energy due to dynamical fric- 
tion against the dark matter background and fall in towards 
the new central object. 

The differential equation for the distance of the satellite 
from the center of the halo (rf r i c ) as a function of time is 
given by 



drfric =-0.428/(6) ^ In A 



( |Binney fc Tremaine 1987| ; |Navarro, Frenk, fc White 1995| ). 
In this expression, m sat is the combined mass of the satel- 
lite's gas, stars, and dark matter halo, and V c is the circu- 
lar velocity of the parent halo. Not to be confused with at 
least two other quantities in this paper denoted by the same 
symbol, here In A is the Coulomb logarithm, which we ap- 
proximate as In A ss ln(l + m^/m^t), where mu is the mass 
of the parent halo. The circularity parameter e is defined as 
the ratio of the angular momentum of the satellite to t hat of 



a circular orb it with the same energy: e = J/J C (E). Lacey 



Cole (1993| ) show that the approximation /(e) 



good approximation for e > 0.02. We draw e for each satel- 
lite from a uniform distribution from 0.02 to 1. A new value 
of e is chosen if the parent halo merges with a larger halo. 

As the satellite falls in, its dark matter halo is tidally 
stripped by the background potential of the parent halo. 
We approximate the tidal radius r t of the satellite halo by 
the condition p Ba t(rt) = Phaio (rfric), i.e. the density at the 
tidal radius equals the density of the background halo at the 
satellite's current radial position within the larger halo. The 
mass of the satellite halo can then be estimated as the mass 
within r t . We assume that both halos can be represented 
by singular isothermal spheres, p oc r~ 2 . When rf r i c is less 
than or equal to the radius of the central galaxy, the satellite 
merges with the central galaxy. 

2.8.2 Satellite- Satellite Mergers 

Satellite galaxies may also collide with each other as they 
orbit within the halo. They may merge or only experience 
a perturbation depending on their relative velocities and in- 
ternal velocity dispersions. From a simple mean free path 
argument, one expects satellites to collide on a time scale 

Tcoll ~ — — (11) 

nav 

where n is the mean density of galaxies, a is the effective 
cross section for a single galaxy, and v is a ch aracteristic ve- 
locity. High-resolution N-body simulations by Makino & Hut 



(1997) indicate that this simple scaling actually describes 



the merger rate quite accurately for collisions of galaxy pairs 
over a broad range of parameter space. They generalize their 
results to obtain an expression for the average time between 
collisions in a halo containing N equal mass galaxies: 



r co n = 500 N' 



( rWoV 



VMpc/ \0.12Mpc 



/ Q-gal \ 4 / (Thalo \ : 

U00kms-V VSOOkms" 1 / 



Gyr. (12) 



dt 



Vc 



(10) 



.lOOkms- 1 / \300kms- 

Here rhaio is the virial radius of the parent halo, r ga i is the 
tidal radius of the dark matter bound to the satellite galaxy, 
tTgai is the internal 1-D velocity dispersion of the galaxy, and 
""halo is the 1-D velocity dispersion of the parent halo. Al- 
though this expression was derived for equal mass galaxies, 
we use it to assign a collision timescale f co ii to each galaxy 
using the mass and tidal radius of each individual sub-halo. 
The probability that a galaxy will merge in a given timestep 
At is then P mr g = At/t C oii- A new velocity dispersion and 
mass is assigned to the post-merger sub-halo by assuming 
that energy is conserved in the collision, and that the merger 
product satisfies the virial relation. Note that we do not al- 
low random collisions between satellite galaxies and central 
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galaxies, even though they are in principle quite likely, be- 
cause we do not know how to model the cross-section for 
such events. 



2.8.3 Merger- Induced Starbursts 

There is considerable observational and theoretical evidence 
that mergers and interactions between galaxies trigger dra- 
matically enhanced star formation episodes known as star- 
bursts. When two galaxies merge according to either of the 
two processes described above, we assume that the cold gas 
is converted to stars at the enhanced rate eburst?n co id/Tdyn, 
where m co id is the combined cold gas of both galaxies, and 
Tdyn is the dynamical time of the larger galaxy. The burst 
efficiency eb ur st ma y depend on the mass ratio of the merg- 
ing galaxies, and is typically between 0.50 to 1. The mean 
properties of galaxies at z = are quite insensitive to the 
details of the treatment of starbursts, although this process 
turns out to be quite important for high redshift galaxies. We 
develop a more detailed treatment of starbursts, ba sed on 



simulations with h y drodynamics and star form ation (Miho: 



fe Hern 



quiot lOOd j- Mihop fe Hornquiot 1006), and invoati 



gate th e implications for high redshift galax ies in a seperate 

paper (somerville, Primack, & Faber 1999). 



2.8.4 Merger-Driven Morphology 

Simulations of collisions between nearly equal mass spiral 



galaxie* produce merger remnants that resemble elliptical 



galaxies. Accretion of a low-mass satellite by a larg er dStc 



will heat and thic ken the disk but not destroy it (Barne; 



fe Herr|quist 19U2 ). However the line dividing these cases 
is fuzzy and depends on many parameters other than the 
mass ratio, such as the initial orbit, the relative inclination, 
whether the rotation is prograde or retrograde, etc. We in- 
troduce a free parameter, /bui ge , which determines whether 
a galaxy merger leads to formation of a bulge component. 
If the mass ratio ms/ms is greater than /bulge, then all 
the stars from both galaxies are put into a "bulge" and the 
disk is destroyed (here ms and ms are the baryonic masses 
of the smaller and bigger galaxy, i.e. the sum of the cold 



gas anc stellar masses). 11 mg/ms < /bui gc , then the stars 
from the smaller galaxy are added to the disk ot the larger 
galaxy. The cold gas reservoirs of both galaxies are com- 
bined. Additional cooling gas may later form a new disk. 
The bulge-to-disk ratio at the observation time can then be 



used to assign rough morphological types. Simien & de Vau 



couleurs (1986) have correlated the Hubble type and the B 



luminosity bulge-to-disk ratio. Using their results, and fol- 
lowing KWG93, we categorize galaxies with B/D > 1.52 as 
ellipticals, 0.68 < B/D < 1.52 as SOs, and B/D < 0.68 as 
spirals. Galaxies with no b ulge are classified as irregulars . 
As shown by KWG93 and [Baugh, Cole, fe Frenk (19964 , 
this approach leads to model galaxies with morphological 
properties that are in good agreement with a variety of ob- 
servations. 



2.9 Stellar Population Synthesis 

Stellar population synthesis models provide the Spectral En- 
ergy Distribution (SED) of a stellar population of a single 



age. These models must assume an Initial Mass Function 
(IMF) for the stars, which dictates the fraction of stars cre- 
ated with a given mass. The model stars are then evolved 
according to theoretical evolutionary tracks for stars of a 
given mass. By keeping track of how many stars of a given 
age are created according to our star-formation recipe, we 
create synthesized spectra for the composite population. A 
free parameter f* um effectively determines the stellar mass- 
to-light ratio; f* um is defined as the ratio of the mass in 
luminous stars to the total stellar mass, mi* um /mt f The 
remainder is assumed to be in the form of brown dwarfs, 
planets, etc. We then convolve the synthesized spectra for 
each galaxy with the filter response functions appropriate 
to a particular set of observations. In this way we obtain 
colours and magnitudes that can be directly compared to 
observations at any redshift. 

Although this approach is satisfying because it results 
in quantities that can be compared directly to observations, 
there are many uncertainties inherent in this component of 
the modelling, as is bound to be the case with such a com- 
plicated problem. The IMF is a major source of uncert ainty. 
The IMF is fairly well determined in our Galaxy ( Scalo 



1986), but we know very little about how universal it is 



or whether it depends on metallicity or other environmental 
effects. The results are somewhat sensitive to the upper and 
lower mass cutoffs as well as the slope of the IMF. Then 
of course there are the difficulties of modelling the com- 
plex physics involved in stellar evolution . Some of the ma- 
jor sources of un certainty mentioned by |Charlot, Worthey" 



Bressan (1996) are opacities, heavy-element mixture, he- 



lium content, convection, diffusion, mass loss, and rotational 



mixing. Compari ng three sets of models, Chariot, Worthey 



& Bressan (1996) find only a 0.05 magnitude dispersion be- 
tween the models in B — V colour, but a larger discrepancy 
of 0.25 magnitudes in V — K colour and a 25% dispersion in 
the mass-to- visual light ratio. However, they also stress that 
there are far greater uncertainties involved in the modelling 
of young (< 1 Gyr) stars, especially stars more massive than 
2M . 

There are currently several versions of stellar population 
models available . We have used the Bruzual fe Chariot (GIS- 
SEL95) models (|Bruzual fe Chariot 1993|; |Charlot, Worthey 



Bressan 1996). These models are for solar metallicity stars 



only. Fo r the results presented in this paper, we have used 
a Scalo (Scalo 1986) IMF and the standard Johnson filters 
provided with GISSEL95. 



2.10 Dust Absorption 

Absorption of galactic light by dust in the interstellar 
medium causes galaxies to appear fainter and redder in the 
ultraviolet to visible part of the spectrum. In this paper, we 
have adopted a sim ple model of dust extinct ion based on the 
empirical results of Wang & Heckman (1996). These authors 
give an expression for the B-band, face-on extinction optical 
depth of a galaxy as a function of its blue luminosity: 



TB 



(13) 



where Ls.i is the intrinsic (unextinguished) blue luminoisity, 
and we use tb,* = 0.8, Lb.* = 6 x W 9 Lq, and (3 — 0.5, as 



© 0000 RAS, MNRAS 000, 000-000 



Semi-Analytic Modelling of Galaxy Formation 9 

Table 2. Cooling/Merging Packages 



name 


cooling 


merging 


starbursts 


Classic 


static halo 


dynamical friction only 


major mergers only; eb urs t = 1 


New 


dynamic halo 


dynamical friction + satellite collisions 


all mergers; e burst = /(mi/m 2 ) 



Table 3. Star Formation/Feedback Packages 



name 


star formation 


feedback 


reheated gas 


Munich 


SFR-M 


eqn. JA ct rn = 2 


all stays in halo 


Durham 


SFR-D 


eqn. pTa r h = 5.5 


all ejected from halo 


Santa Cruz (fiducial) 


SFR-M 


disk/halo 


disk/halo 


Santa Cruz (high fb) 


SFR-M 


disk/halo 


disk/halo 


Santa Cruz (C) 


SFR-C 


disk/halo 


disk/halo 



found by Wang & Heckman (1996). We then relate the B- 
band optical dept h to other bands using a standard Ga lactic 
extinction curve jCardelli, Clayton, fc Mathis 1989| ). The 
extinction in magnitudes is then related to the inclination 
of the galaxy using a standard "slab" model (a thin disk 
with stars and dust uniformly mixed together): 



-2.5 log 



— T\ sec 8 



t\ sec 9 



(14) 



wher e 6 is the angle of inclination to t he line of 



sight (Guiderdoni fc Rocca-Volmerange 1987; Fully & 



Fouquc 1985 ) . We assign a random inclination to each model 
galaxy. The extinction correction is only applied to the disk 
component of the model galaxies (i.e. we assume that the 
bulge component is not affected by dust). 

2.11 Model Packages 

There are many possible permutations of the different in- 
gredients that we have introduced above. In the interest of 
practicality, we have chosen several "packages" of ingredi- 
ents to explore in this paper. In the relevant sections, we 
comment on which elements of the package are important 
in determining various quantities. 

In order to understand the effects of the new ingredients 
that we have introduced pertaining to cooling, galaxy merg- 
ing, and starbursts, we introduce two "cooling/merging" 
packages (see Table ^). The ingredients of the "Classic" 
package were chosen to be similar to the published Mu- 
nich and Durham models In this paper, we always ap- 
ply them within the SCDM cosmology, which was used for 
much of the previous work. The ingredients of the "New" 
cooling/merging package reflect additions or modifications 
to our models. We apply the "New" models within more re- 
alistic (or anyway more fashionable) cosmologies. The new 
ingredients are describ ed a bove, namely, the dynamic -hal o 
cooling model (Section 2.3), satellite collisions (Section * 



and more detailed modelling of starbursts in galaxy-galaxy 

T Except that the cooling model used in the published Durham 
mod els is more similar to our "dynamic halo" model (see Sec- 
tion [?.3| ). 



mergers. In the "Classic" models, starbursts occur only in 
major mergers and with efficiency eb urs t = 1- In the "New" 
models, starbursts occur in all mergers and eburst is a func- 
tion of the mass ratio of the merging galaxies. The details of 
the starburst modelling are of minor importance for galaxy 
properties at z ~ 0, and will be dealt with in detail in a 



Faber 1999) 



companion paper (Somerville, Primack, 

We also wish to understand the effects of different 
choices of star formation and supernova feedback recipes, 
and introduce several "sf/fb" packages. We have chosen the 
ingredients of the first two packages to be similar to the 
choices made by the Munich and Durham groups with re- 
spect to the star formation and supernov a fee dback, includ- 
ing the fate of reheated gas (see Section 2.6). It should be 
kept in mind, however, that although we will refer to these 
as the "Munich" and "Durham" packages, our models differ 
from those of these other groups in many respects and we 
are not trying to reproduce their results in detail. On the 
contrary, we wish to isolate the effects of the way that star 
formation and supernova feedback are modelled. For exam- 
ple, as we will discuss in Section |^, the published models of 
the Munich and Durham groups are normalized such that 
a galaxy of a given circular velocity is considerably fainter 
than in our models. Here we will always normalize all of the 
packages in the same way, as described in Section ^. We will 
refer to the results that we obtain from our code, normal- 
ized as described in this paper, as the Munich and Durham 
"packages". When we wish to refer to the results obtained 
by the Munich and Durham groups using their codes, nor- 
malized in their own ways, we shall refer to the "actual" or 
"published" Munich or Durham models. 

The third package, which we refer to as the "Santa Cruz 
(fiducial)" package, is a hybrid of the Munich-style star for- 
mation law (SFR- M) a nd the disk-halo feedback model de- 
scribed in Section |2.5| . "Santa Cruz (high fb)" is the same 
as Santa Cruz (fiducial) except that the supernova feedback 
parameter is turned up by a factor of five. The "Santa Cruz 
(C)" package assumes that the star formation efficiency is 
constant at all redshifts and in galaxies of all sizes (SFR-C). 
Note that this is equivalent to the star formation law sug- 
gested for use with milder supernova feedback (/ v = 0.01) 
by CAFNZ94. Again, we combine this with the disk-halo 
feedback model. 
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Zhot (2.3): the assumed metallicity of the hot gas, used 
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Figure 4. Fits to the observed I-band Tully-Fisher relation from 
several samples. PT92 is the (Pierce & Tully 1992) relation for 
nearby galaxies, scaled assuming variou s values of the Hubble 



Giovanelli et al. (1997) 



parameter. Giovanelli et al. is from the 
sample of cluster spirals. HM and MAT are the Ha n-Mould an d 
Mathcwson et al. samples discussed in Willick et al. (1995, 199C). 
The magnitude of a "Milky Way" galaxy (V c = 220 km s _1 ) in 
the Durham models, from Fig. 11 of CAFNZ94 is indicated by 
the open square. The approximate magnitude of a "Milky Way" 
galaxy in the Munich models is indicated by the open triangle. 
The magnitude of the "Milky Way" in our models is indicated by 
the filled hexagon. 



in the cooling function. 



3.2 Free Parameters 

The free parameters and the sections in which they were 
introduced are: 

the star formation timescale 
supernova reheating efficiency 
y (12. 7f) : chemical evolution yield (mass of metals pro- 
duced per uni t mass of stars) 

f* um (2.9): the fraction of the total stellar mass in lu- 




minous stars 

• /mrg 



(2.8.1) 



the initial distance of satellite halos from 
the central galaxy after a halo merger, in units of the (post- 
merger) virial radius. 

• /bulge (sec:models:morph): the mass ratio that divides 
major mergers from minor mergers; determines whether a 
bulge component is formed 

To set the values of the free parameters, we define a 
fiducial "reference galaxy", which is the central galaxy in a 
halo with a circular velocity of V c = 220 km s~ . We set the 
most important free parameters by requiring the properties 
of this reference galaxy to agree with observations for an 
average galaxy with this circular velocity. As an important 
constraint, we would like to require our reference galaxy to 
have an average I magnitude given by the Tully-Fisher re- 
lationship (we normalize in I rather than B because it is 
less sensitive to recent starbursts and the effects of dust). 
But first, we discuss a subtlety in the process of compar- 
ing the models with this observation, which has led to some 
confusion in the past. 



3 SETTING THE GALAXY FORMATION 
PARAMETERS 

We have introduced a number of parameters. "Free" param- 
eters are adjusted for each choice of cosmology and package 
of astrophysical recipes, according to criteria that we shall 
describe. The "fixed" parameters keep the same values for 
all of the models. In this section, we summarize the param- 
eters and the procedure that we use to set them. 



3.1 Fixed parameters 

The physical meaning and values of the fixed parameters, as 
well as the section in which they are discussed in detail, are 
summarized as follows: 



3.2.1 Tully-Fisher Normalization 



a a* = —1.5 (2.5): power used in power-law (Durham) 
star form ation law 

• a r h (2.6.1): power used in supernova reheating power- 
law (a r h = 2 for Muni ch, a r h = 5.5 for Durham) 

• /rh.disk = 2 ( |2.6,2 ) : fudge factor used in disk-halo feed- 
back model 



• /reheat = 0.5 fl2.3.2| ): In the dynamic halo cooling model, 
the cooling time r coo i is reset after a halo merger event if the 
largest progenitor of the current halo is less than a fraction 

/reheat of its HiaSS. 



If we use a local sample, such as that of fierce fc Tully 
(1992| ) ; the relation between absolute magnitude and line- 



width has been determined by measuring a distance to each 
galaxy using various standard methods (e.g. Cepheids, RR 
Lyrae, planetary nebulae). This relation therefore i ntrinsi- 



cally contains an effective Hubble parameter. The Pierce 
k. Tully (1992 ) sample, when used to derive distances to 
the Ursa Major a nd Virgo clusters, im plies Ho = 85 ± 
10km s" 1 Mpc" 1 ( [Pierce fc Tully 1988; ). Our models are 
set within predetermined cosmologies with various values of 
the Hubble parameter (H = 50, 60, or 70 km s" 1 Mpc -1 ). 
To compare this data to the different cosmologies used in 
our models, one approach is to simply scale the observed 
absolute magnitudes: 



Afmodol — Mobs + 5 log(/l mo del/^obs ) 



(15) 



This is effectively what CAFNZ94 say they have done in 
their Fig. 11, assuming h OOB — 1.0 (although it looks more 
as though they used h = 0.80). They then interpret their 
Fig. 11 as indicating that their models are discrepant with 
the observed Tully-Fisher relation because their galaxies are 
~ 1.8 magnitudes too faint at a given circular velocity. In 
contrast, KWG93 claim good agreement with the Tully- 
Fisher relation, and show this in their Fig. 7. This leaves 
one with the impression that a Durham galaxy would be 
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Table 4. Galaxy formation parameters for the "Classic" models 
(SCDM) 



model 


i~2 


e° 
SN 


y 


f* 

J lum 


/bar 


Munich 


100 


0.2 


1.3 


1.0 


0.1 


Durham 


4.0 


N/A 


1.8 


1.0 


0.125 


Santa Cruz (fiducial) 


100 


0.1 


1.8 


1.0 


0.125 


Santa Cruz (high fb) 


100 


0.5 


4.0 


1.0 


0.125 


Santa Cruz (C) 


8.0 


0.125 


1.8 


1.0 


0.125 



Table 5. Galaxy formation parameters for the "New" (Santa 
Cruz fiducial) models 



model 




e° 
e SN 


V 


f 

J lum 


/bar 


SCDM 


100 


0.1 


1.7 


1.0 


0.125 


tCDM 


100 


0.05 


l.S 


1.0 


0.11 


ACDM.5 


100 


0.125 


3.5 


1.0 


0.11 


OCDM.5 


100 


0.125 


3.5 


1.0 


0.11 


ACDM.3 


50 


0.125 


2.2 


0.80 


0.13 


OCDM.3 


80 


0.125 


3.7 


0.9 


0.13 



about 1.8 magnitudes fainter than a Munich galaxy with 
the same circular velocity. 

It is difficult to make a direct comparison from the 
published papers because KWG93 plot the TF relation in 
the B band, and in terms of luminosity, whereas CAFNZ94 
use the I band, and plot Mi — 5 log h. However, we can 
easily check what would happen if we scaled the B-band 
relation used by KWG93 in the same way. If we assume 
^B.obs = —20.7 for a V c = 220 km s _1 galaxy, from the 



Pierce & Tully (1992) relation, and take h ^ s = 1.0 and 
^modci = 0.5, this would imply Mr^modci = —22.2. This is 1.2 
to 2.2 magnitudes brighter than the "Milky Way" normal- 
ization (Mb ~ —20 to -21), used in the published Munich 
models. What this means is that the apparent good agree- 
ment with the TFR seen in Fig. 7 of KWG93 is because they 
assumed h obs = h mo dei = 0.5. 

One lesson in all of this is that trying to normalize the 
models with the local Tully-Fisher data is problematic be- 
cause it is so sensitive to the Hubble parameter. A more 
robust approach is to use the velocity-based zero-point from 
the compilati on of sever al mo re distant TF sample s from 
Willick et al. ( |l995| , |l996|) and |Giovanelli et al. (19971 ). This 
effectively gives us a relation between Alj — 5 log h and line- 
width. The Hubble parameter is explicitly scaled out, so we 
can apply the normalization fairly across models with dif- 
ferent values of Ho- We show the fits to these observed rela- 
tions along with the I magnitude of a "Milky Way" galaxy 
for the published Munich and Durham models and for our 
fiducial models in Fig. |^. The magnitude of the Durham 
"Milky Way" is taken from Fig. 11 of CAFNZ94. The Mu- 
nich "Milky Way" is the I-band magnitude that we get in 
the version of our code in which we try to reproduce all the 
assumptions of the Munich models, and set the free param- 
eters to get the B magnitude they quote (Mb ~ —20). It is 
approximately the same as in the actual Munich models (G. 
Kauffmann, private communication). We can now see that if 
placed side by side, the Munich and Durham model galaxies 
have almost the same magnitudes at V c = 220 km s" 1 , and 
that both are about 2 magnitudes fainter than the observed 
I-band TFR, independent of assumptions about the Hubble 
parameter. This can be reconciled with the Munich group's 
assertion that they reproduce the observed B-band TFR by 
two factors. One is the scaling with Ho that has already been 
discussed. The second factor is that the mo del galaxies are 
too blue in B-I. We also see that the local Pierce fc Tully 
(1992Qelation only agrees with the results of more distant 
samples if a relatively high value of the Hubble parameter 
(h ~ 0.85) is assumed. This is further evidence that some 



sort of rescaling is necessary if the local relation is to be 
used in conjunction with theoretical models. 

Therefore the results of the published Munich and 
Durham models are actually more consistent than it ap- 
peared. The reason their "Milky Way" is so much fainter 
is easy to understand. The value of the parameter that we 
refer to as f* um is 0.5 in the Munich models and 0.37 in the 
Durham models. This corresponds to the assumption that 
50% or 63% of the stellar mass is in the form of non-luminous 
brown dwarfs or planets. This is an unrealistically large con- 
tribution from non-luminous stars according to most theo- 
ries of star formation, and results in stellar mass-to-light 
r atios about a factor of 2-3 higher than the observed values 
( |Gilmore 1997 ). Taking f* uln — 1 results in more reasonable 



stellar mass-to-light ratios, and brings the reference galaxy 
into better agreement with the TFR. 



3.2.2 Setting the Free Parameters 

We now set our parameters to get the central galaxy in a 
V c = 220 km s^ 1 halo to have M/- 5 log ft ~ -21.6 to -22.1, 
which is consistent with the values predicted from the fit- 
ted relations for the three distant I-band surveys mentioned 
above. Because the observations have been corrected for dust 
extinction, we use the non dust-corrected magnitude of the 
reference galaxy to normalize the models (actually we should 
use magnitudes with face-on dust corrections, but in the I- 
band these are quite small). To convert between the mea- 
sured Hi line-widths Wr and the model circular velocities, 
we assume Wr = 2V C . However, it should be kept in mind 
that this transformation is not necessarily so straightfor- 
ward, and this could change the slope and curvature of the 
relation especially on the small- linewidth/faint end. Note 
that we have also implicitly assumed that the rotation ve- 
locity of the galaxy is the same as that of the halo, i.e. that 
the rotation curve is flat all the way out to the virial ra- 
dius of the halo. This neglects the effect of the concentrated 
baryons in the exponential disk, which will increase the cir- 
cular velocity at small radii. Moreover, if the dark matter 
halo profiles resemble the form found by Navarro, Frenk 



8z White (1996), for galaxy-sized halos, the rotation veloc- 



ity at ~ 2 disk scale lengths (where the observed TFR is 
measured) is 20-30% higher than that at the virial radius of 
the halo. This would mean that our V c = 220 km s" 1 galaxy 
would live inside a V c ~ 180 km s _1 halo. Both of these ef- 
fects would lead to a larger galaxy circular velocity for a 
given halo mass, hence to smaller mass-to-light ratios. 

We also require our average reference galaxy to have a 
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Figure 5. The history of stars, cold and hot gas within all ha- 
los that will eventually form a "local group" (V c = 220km s — 1 ) 
sized halo at z = in the Classic/SCDM models. The solid lines 
indicate the total baryon fraction (stars + cold + hot) within the 
halo, with respect to the universal value. Dotted, short, and long 
dashed lines indicate the fraction of baryons in the form of stars, 
cold gas, and hot gas respectively. 



cold gas mass of m coW ~ 3.2 x 10 9 to 2.5 x W 10 h~ 2 M Q . 
This is consistent with the average Hi mass of a galaxy with 
Mr - 5 log h ~ -21.8 (|de Blok, McGaugh, fe van der Hulst 



1996|)~fcul tiplied by a factor of two to account for molecular 
hydrogen ( [Young fc Knezek 1989 ). This fixes the two main 



parameters r" and e'g N , although there is some unavoidable 
degeneracy (see Section ^). The yield y is set by requiring 
the stellar metallicity of the reference galaxy to be equal 
to solar. Note that the value of y does not affect any other 
properties of the galaxies b ecau se we have used a fixed hot 
gas metallicity (see Section 2.3 ). 

Following Baugh, Cole, & Frenk (1996a), the parame- 
ter /bulge is fixed by requiring the fraction of morphological 
types to be approximately E/S0/S+Irr = 13/20/67 (these 
ratios were obtained by scaling the observations of Love- 
day (199(i) to account for unclassifiable galaxies). Using 
/bulge = 0.25 results in roughly these fractions for all the 
models investigated here, and this value is used throughout 
this paper. 

The values of the free parameters used in the models 
presented in this paper are given in Table ^ and Table ^[ 
We run many realizations and use the average values of these 
quantities in order to fix the values of the free parameters. 



4 THE FORMATION OF AN L, GALAXY 

As we have discussed, we normalize our models by requiring 
certain properties of a reference galaxy about the size and 
luminosity of the Milky Way to agree with observations. In 
this section, we illustrate how the formation history of our 
reference galaxy and its satellite companions depends on 
the prescriptions we use for star formation and supernova 
feedback (sf/fb), and the values of our free parameters. This 
will help in interpreting the results of the next section, in 
which we show how global quantities such as the luminosity 
function depend on these assumptions. 

Table |i| shows the fiducial values of the free parame- 
ters used for each of the sf/fb packages introduced in Sec- 
"Classic" cooling/merging package and the 



tion 2.11 



The 



SCDM cosmology are used for all of these models, and the 
dynamical friction parameter is set to /mrg = L Table || 
shows the parameters used for the models with the "New" 
cooling/merging package and the Santa Cruz (fiducial) sf/fb 
package. For these models, we set / mrg = 0.5, which is in 
better agreement with the results of high resolution N-body 
simulations (A. Klypin, private communication). It should 
be noted that the specific values of these parameters may de- 
pend somewhat on the details of the implementation of our 
code. Also note that r* and eg N do not function in precisely 
the same way in the different packages because of the dif- 
fering functional forms of the recipes, so they cannot always 
be compared directly. 

Fig. |B] illustrates the redshift evolution of the baryonic 
content (stars, cold gas, and hot gas) of halos that will even- 
tually form a "local group" (V c = 220 km s _1 ) sized halo at 
z = 0. The dependence on the sf/fb package and the value of 
the supernova feedback parameter is shown. Note that star 
formation occurs much earlier in the Munich package than 
in the Durham package models. This is due to two combined 
effects. First, as we discussed in Section 2.5, with SFR-M the 



star formation efficiency is higher at high redshift because 
the typical galaxy dynamical times are shorter. In SFR-D, 
star formation is less efficient in objects with smaller cir- 
cular velocities. At high redshift the characteristic circular 
velocities tend to be smaller, so this leads to less star forma- 
tion. Second, the much stronger supernova feedback in the 
Durham package models leads to additional suppression of 
star formation, especially in small objects. 

The bottom four panels break down these ingredients. 
In the Santa Cruz (fiducial) package, we use SFR-M. The 
effect of turning up the feedback efficiency by a factor of five 
is shown in the right panel. Star formation is suppressed, 
and more so at higher redshift where objects are smaller, 
but the effect is not as dramatic as in the Durham package. 
The bottom-most panels show the Santa Cruz (C) package, 
which assumes constant star formation efficiency (SFR-C). 
This package is intermediate between the Durham package 
and the Santa Cruz SFR-M (fiducial) package. 

Fig. ^ shows how tuning the free parameters changes 
the properties of the reference galaxy in the "Classic" Santa 
Cruz (fiducial and C) models, within the SCDM cosmology. 
The figure shows the space of I-band magnitude and cold 
gas mass, along with the target area used to normalize the 
models (shaded box). Symbols show the location of average 
reference galaxies within this space for different values of the 
free parameters t*, eg N , and /baryon - The dependence on the 
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Figure 6. The change in I-band magnitude and cold gas mass 
of an average reference galaxy in the Classic/Santa Cruz models 
(top panel: fiducial; bottom panel: SFR-C), SCDM cosmology, 
as the free parameters are tuned. As is increased, one moves 
upwards along the solid lines connecting symbols of the same 
shape. As 6g N is increased, one moves leftwards along the dotted 
lines and the symbol shapes change. This "grid" is run with fixed 
/baryon = 0.076. The dotted line connecting filled triangles (top 
panel) shows the effects of varying the baryon fraction /baryon i 
with fixed r° = 100 and e gN = 0.25. 



Figure 7. The change in I-band magnitude and cold gas mass 
of an average reference galaxy in the New/Santa Cruz (fiducial) 
models, rCDM and ACDM.3 cosmologies, for different values of 
the free parameters. Open symbols connected by dotted and solid 
lines represent the same values of t® and eg N as in Fig. ^, with 
/baryon = 0.076,0.129 (rCDM, ACDM.3) and /* um = 1. Filled 
triangles in the top panel represent varying the baryon fraction 
as in Fig. H (top panel). In the bottom panel, the filled triangles 
represent varying values of = 0.2, 0.4, 0.6, 0.8, 1.0 (left to 

right). 



free parameters takes a different form for different star for- 
mation/feedback recipes. Generally, increasing the star for- 
mation timescale r* leads to an increased gas mass, and to 
a much lesser extent, a fainter luminosity. Increasing the su- 
pernova feedback efficiency eg N leads to a fainter luminosity 
and, to a lesser extent, smaller gas mass. Increasing /baryon 



leads to a larger gas mass and luminosity. The same exer- 
cise is repeated in Fig. [?] with the New/Santa Cruz (fiducial) 
package, for the rCDM and ACDM.3 cosmologies. Here we 
show the effect of varying f* um , which can only make galax- 
ies fainter as it can only take values less than one. Figure ^ 
shows the location of the average reference galaxy in this 
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Figure 8. Average I-band magnitude and cold gas mass of a 
fiducial reference galaxy in the Santa Cruz models. Error bars 
indicate the standard deviation in these quantities over many 
realizations. 



space and the standard deviation of these quantities over 
many ensembles, for all the Santa Cruz models, for the final 
fiducial values of the free parameters shown in Tables ^| and 



We would have liked to consider /baryon to be de- 
termined independently, thus eliminating a free parame- 
ter. However, it is apparent from Fig. |6| that if we take 
/baryon = 0.076, which corresponds to the baryon fraction 
derived from observations of de uterium at high redshift, 
fl b h 2 = 0.019 (fTytler et al. 1999|) for h = 0.5 and fi = 1, 



the reference galaxy is too faint and gas poor in the fio = 1 
cosmologies compared to our desired normalization. The val- 
ues of /baryon — 0.11 to 0.13 that we find necessary to obtain 
our desired normalization are similar to those typically de- 
rived from_j^rv_diffe^ in groups and clus- 
ters ( Mohr, Ma thicscn, fc Evrard 1999 ). As emphasized by 
White et al. (1993), for high values of Qp ~ 1 this is i n- 



consistent with Big Bang Nucleosynthesis ( Copi, Schramm 
& Turner 1996), and with the measurement of Tytler et al 
(1999^7 This could be interpreted as further evidence from 
the galaxy side that f2o is probably less than unity. It is cu- 
rious that the best agreement occurs for Qo — 0.4 — 0.5 and 
h = 0.6 — 0.65, just the values currently favored by indepen- 
dent considerations. But given the large and numerous un- 
certainties in our modelling (particularly cooling, feedback, 
star formation efficiency, and the IMF), we do not regard 
this as much more than a curiousity, albeit a rather com- 
fortable one. For example, we have neglected the eventual 
return of the gas expelled by supernovae, and the recycled 
gas from dead stars. If these were included, we might be able 
to reduce the value of /baryon somewhat. For the moment, 
we formally consider /baryon and f* nm to be simply free pa- 
rameters, which are close enough to their plausible physical 
values as to not cause too much concern. 



Figure 9. The mass function of dark matter halos predicted by 
the standard Press-Schechter model for various CDM cosmolo- 
gies (light broken lines). The bold lines show the mass function 
of galactic halos, estimated from the observed APM luminosity 
function as described in the text, for an SCDM or tCDM cosmol- 
ogy (dotted), and for the ACDM.3 cosmology (long dashed-dotted 
line); other cosmologies lie between these two cases). The short 
solid line shows a power-law with slope a = —2. 



5 COMPARISON WITH LOCAL 
OBSERVATIONS 

In this section, we investigate the predictions of our models 
for a number of important galaxy properties. We have sev- 
eral goals: we compare the results of our models with previ- 
ously published work, explore the importance of the choice 
of sf/fb recipe and the values of certain free parameters, and 
compare with observational results. 

The local number density of galaxies as a function of 
their luminosity is clearly a key prediction of any success- 
ful model of galaxy formation. The halo mass function pre- 
dicted by any of the currently popular CDM-based models 
has a very different shape from the characteristic Schechter 
form of observed luminosity functions. At masses less than 
10 13 h- 1 M , the CDM mass function is a power law with 
a slope a ~ — 2, much steeper than the faint-end slope of 
the observed field galaxy luminosity function a ~ —1.0 to 
— 1.5. The exponential cut-off occurs at ~ 10 14 h" 1 Mo, 
much larger than the expected halo mass corresponding to 
an L* galaxy. In Fig. |^, we show the halo mass function pre- 
dicted by the standard Press-Schechter model, along with 
the mass function of galactic halos estimated in the follow- 
ing simple way. We find the circular velocity of the dark 
matter halo associated with a typical L* galaxy using the 
Tully-Fisher relation and assumin g the rotation curve i s flat 
19.5 (cf. ~ 



Loveday et al. 1992 ) , and 



: Tully 1992; Tully 



Using Mg — 5 log h = 
the B-band Tully-Fisher relation ( Pierce 
et al. 1997) scaled to h ^ B 

V* = 160 km s" 1 . We can then translate this to a mass 
using the spherical top-hat model (see Appendix A). From 



0.80 (see Section fl), we find 
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Figure 10. The B-band luminosity function of galaxies for the 
Classic/SCDM models. Crosses and open squares show the mod- 
els with the original Press-Sc hechter weighting and t he improved 
Press-Schechter weighting of riheth & Tormcn (1999), both with- 
out dust extinction. Solid squares show the models with inclusion 
of the empirical dust models and the improved Press-Schechter 
model. Dashed lines indicate the fits to the observed luminos- 
ity function from several redshift surveys as indicated in the key 
(references given in the text). 



Fig. 



Al 



we see that this corresponds to a halo with a mass 
of about 1.0 x 10 12 M Q to 1.3 x 1O 12 M , depending on the 
cosmology. Using this constant light-to-mass conve rsion, wc 
transla te th e observed B-band luminosity function (Loveday 
et al. 11)92) to the mass function shown in Fig. ^| Of course 
this translation is complicated by sub-structure (each of the 
halos in the Press-Schechter model may contain multiple 
galaxies of various sizes) , as well as by the varying mass-to- 
light ratio of galaxies of different morphological types and 
other complications. However, these effects will introduce 
changes of order a factor of a few, and the discrepancy is 
much larger. The problem may be summarized as follows. 
In order to get from any CDM mass function to the observed 
luminosity function, it seems that the conversion from halo 
mass to galaxy luminosity must be more complicated than 
what we have assumed in this simple calculation; in partic- 
ular, apparently the mean mass-to-light ratio must decrease 
as we move away from V* in both directions. On the other 
hand, the constant mass-to-light model gives us a perfect 
power-law Tully-Fisher relation with the correct slope and 
zero scatter. Any scatter in the mass-to-light ratio at fixed 
V c will introduce scatter in the TFR, and any systematic 
variation with V c will introduce curvature. Satisfying both 
constraints simultaneously has proven to be a challenge. 
For example, the first generation of Munich and 



-1 
-2 

-3 
-4 



7 -1 
o 

CO 

£ -3 
o _ 4 

W) 


-1 



_ 1 1 1 1 1 1 

r SCDM 

" 


iiii 


_ 1 1 1 1 1 1 1 1 1 1 1 1 j B l !_ 

: tCDM .. 

:: / ~ 
It ■ 

r T ' Z 


r ACDM.5 

r / 

- a t 

□ m 

□ 1 


■ "J 


: 0CDM.5 t , ■ ■ * 


7 ACDM.3 

: j 

- fm 

□ B 

- n □ ! ^ 


■ 

1 1 1 1 1 1 1 1 1 - 


7 0CDM.3 

A 

□ M 

~ n m _Z 
□ □ ■// 



-22 -20 -18 -16 



-22 -20 -18 -16 



M B - 5 log h 



Figure 11. The B-band luminosity function of galaxies for the 
New/Santa Cruz (fiducial) models. Key as in Fig. [Io| 



Durham models effectively assumed a stellar mass-to-light 
ratio a factor of 2-3 times larger than the face-value pre- 
diction of the Bruzual- Chariot models (i.e. f{ um — 0.5 in 
the Munich models and 0.63 in the fiducial Durham mod- 
els). This pushed the galaxy mass function (bold curves in 
Fig. to the right, to the point where the number den- 
sity roughly agreed at the "knee" (L*). However, it made 
the galaxies about 2 magnitudes too faint compared to the 
observed Tully-Fisher relation. The Durham group designed 
their star formation and supernova feedback models in order 
to obtain light-to-mass ratios that decreased rapidly with 
V c - This flattened the faint end slope of the luminosity func- 
tion but led to a pronounced deviation from the observed 
power- law shape of the Tully-Fisher relation (cf. Fig. 11 of 
CAFNZ). In the following two sub-sections we discuss our 
results for these two fundamental observed quantities. 

5.1 The Luminosity Function 

We show the B-band luminosity functions for the Classic 
SCDM models in Fig. [uj (the packages are summarized in 
Table . The curves show fits to the observed B-band lumi- 
nosity funct ions deri ved from the CfA (Marzke, Huchra, & 
Gcllcr 19941 ), APM ( |Loveday et al. 1992| ), SSRS ( |da Co~ 
et al. 1994 ), ESP ( jZucca et al. 1997|), UKST ( |Ratcliffc 
et al. 1997| ), and 2dF ( |Folkcs et al. 1999] ) redshift surveys. 
The observational fits have been converted to the John- 
son B filter band used in our models using t he conversion 
Mb T = Mz — 0.45 for Zwicky magnitudes ( Shanks et al 
M„ 



1984 ) and m b = M bj + 0.2 (CAFNZ94). We show the ef- 
fects of using the improved Press-Sch echter weighting from 
the model of Sheth fc Tormen (1999 ), and of correcting for 
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Figure 12. The K-band luminosity function of galaxies for the 
Classic/SCDM models. Crosses and open squares show the mod- 
els with the original Press-Sc hechter weighting and t he improved 
Press-Schechter weighting of rlheth & Tormcn (1999), both with- 
out dust extinction. Solid squares show the models with inclusion 
of the empirical dust models and the improved Press-Schechter 
model. Dashed lines indicate the fits to the observed luminosity 
function as indicated in the key (references given in the text). 



dust extinction using the recipe discussed in Section 2.10 
Clearly both of these effects help to alleviate the tendency 



of the r lodels to overpredict the number density of galaxies. 
The extinction correction is larger for luminous galaxies (as 
a direct result of the Wang-Heckman recipe), but recall that 
the correction is only applied to the disk component of our 
galaxies. Early type galaxies (which are defined as having 
large bulge-to-disk ratios) therefore suffer much smaller cor- 
rections. It appears plausible that extinction due to dust is 
an important factor in reconciling the discrepancy between 
the modelled B-band luminosity function and Tully-Fisher 
relation. It should be noted that the observed B-band lumi- 
nosity function derived from any of the above redshift sur- 
veys is not corrected for the effects of dust extinction, and 
Tully-Fisher work always includes a correction for both in- 
ternal and Galactic dust extinction. This has been ignored in 
the previous theoretical comparisons that we have discussed. 
As a point of reference, note that the correction for i nternal 
dust ex tinctio n in M31 ranges fro m 0.27 magnitudes ( Pierce 
& Tull"~ 



1992[ ) to 1.0 magnitude ( |Bernstein et al. 1994] ) 



the I-band. This is to stress that both the corrections and 
the uncertainties associated with dust extinction are large. 
The corrections are presumably even larger in the B-band, 
and for more inclined galaxies. 

However, in our models, dust extinction has a negligible 
effect in faint galaxies, and very strong feedback (Durham or 
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Figure 13. The K-band luminosity function of galaxies for the 
New/Santa Cruz (fiducial) models for different cosmologies. Key 
as in Fig. [l^ , 

Santa Cruz with high feedback) still seems to be necessary to 
reproduce the observed faint-end slope within SCDM. The 
Santa Cruz SCDM models with more moderate feedback 
produce a factor of ~ 3 — 4 excess of SMC-sized galaxies, 
which is probably difficult to reconcile with observations, 
even accounting for sources of incompleteness such as surface 
brightness selection effects ( Dalcanton, Spergel, & Summers 



1997] ; palcanton et al. 19971 ). 

Fig. |ll] shows the luminosity function of the other cos- 
mological models, using the New/Santa Cruz (fiducial) pack- 
age. The tCDM model looks quite similar to the SCDM case 
and similarly shows an excess of faint galaxies. We have tried 
several variations of the rCDM model shown here in an at- 
tempt to correct this. One might think that lowering the 
normalization as would decrease the overall number density 
of galaxies. Actually, as mainly controls the location of the 
exponential cut-off in the mass function. As we showed in 
Fig. ^j, this lies well above the scale of galactic halos, and so 
changing as within the bounds allowed by the observed clus- 
ter abundance does not significantly improve our results. We 
also tried to reduce the number of faint galaxies by increas- 
ing the merging rate (we decreased the dynamical friction 
merging timescale to /mr g = 0.1), but we find that this does 
not improve the faint-end significantly and leads to a severe 
excess of bright galaxies. 

However, the low-fi models, particularly the fin = 0.5 
models, reproduce the overall shape and normalization of 
the observed luminosity function remarkably well. Note that 
at very faint magnitudes (Mb — 51og/i > —17) the ob- 
served luminosity function is not well determined, but there 
is actually a suggestion of the steepening faintwards of 
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Figure 14. The Tully-Fisher relation for the Classic/SCDM mod- 
els. Broken lines show fits to the observed I-band Tullv-Fisher re - 
lation from several samples (G97 is from |5iovanelli et al. (1997] ). 
HM is the Han-Mould s ample, and MAT is Mathewson et al. sam- 
ple from Willick et al. (1995, 1996 )). The bold horizontal dashed 
line shows the approximate magnitude limit of the observations. 
The symbols show the results of the models (crosses show cen- 
tral galaxies only and filled squares show central and satellite 
galaxies), and the error bars indicate 1-cr variances over different 
merger history realizations. Only model galaxies that contain cold 
gas and are identified as spirals are included. 



Mb — 51ogfo > — 1 7 that we see in our Qp = 0.5 models 



Marzke, Huchra, & Geller 1994 



Folke; 



( Zucca et al. 1997 
et al. 1999). The fit on the bright end could be improved 



by adjusting the parameters of our dust recipe, which w e 



have t aken at face value from Wang fa Hockman (1006 ^- 
The f2pl— 3 models show a slight deficit nf galavies around 
L*, even without dust extinction, but this is within the un- 
certainties on the normalization of the observed luminosity 
function. 

The effects of dust are significantly reduced in longer 
wavelength bands such as the near IR, however the observed 
luminosity function is not as well determined as it is in op- 
tical bands. We compare our results with two recent deter- 
minations of the K-band luminosity function (we use the K s 
filter (referred to as K' in the IRIM manual) downloaded 



from the KPNO website, ftp://ftp.noao.edu/kpno/filterf, 
with standard Vega zeropoints). The Classic/SCDM models 
are shown in Fig. I14 The wide-field K-band survey discussed 



in Gardner et al. 7| »' »7| ) covers an area of ~ 4.4 square de- 



grees, and probably provides the best existing determination 
of the bright end o f the K-band lumin osity function. The 
survey discussed in Szokoly et al. (199S| ) has a smaller area 
(0.6 square degrees) but has a fainter limiting magnitude, 



Figure 15. The Tully-Fisher relation for the New/Santa Cruz 
(fiducial) models. Broken lines show fits to the observed I-band 
Tully-Fisher relation from several samples (see Fig. |wj). The bold 
horizontal dashed line shows the approximate magnitude limit 
of the observations. The symbols show the results of the models 
(crosses show central galaxies only and filled squares show central 
and satellite galaxies), and the error bars indicate \-a variances 
over different merger history realizations. Only model galaxies 
that contain cold gas and are identified as spirals are included. 



and thus presumably provides a more reliable estimate of 
the faint-end slope. 

All the SCDM models show an overall excess of galax- 
ies of all luminosities, and the Munich package shows a 
slightly steeper faint-end slope than the observations. Both 
Santa Cruz packages and the Durham packag e have a faint 



end sl ope consistent with the observations of szokoly et al 



(1999). The New/Santa Cruz (fiducial) models for the other 
cosmologies are shown in Fig. 13 The fiducial rCDM mod- 
els now show a good match on the bright end but still have 
an excess on the faint end. The flo = 0.5 models are a near 
perfect fit over the range of luminosities probed by the ob- 
servations, except in the very brightest bins. They do not 
cut off as sharply as a pure Schechter function at brighter 
luminosities, but the observed luminosity function is not well 
determined on the bright end because of small samples and 
evolutionary effects. It should be noted that the evolutionary 
and k-corrections applied to the data are non-negligible, and 
are cosmology dependent. The fits shown here are for a Uni- 
verse with qo = 0.5, which is inconsistent with our low-fio 
cosmologies. A more detailed comparison with the observa- 
tions is clearly in order; however, given these uncertainties 
the level of agreement shown here is encouraging. 
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5.2 The Tully-Fisher Relation 

Recall that we have adjusted the free parameters in our 
models to force our reference galaxy to lie on the I band 
Tull y-Fisher relation deriv ed by Willick et al. ( J.994 |l99fj ) 
and jGiovanelli et al. (1997| ). Fig. shows the fits from the 
three observational samples mentioned above and the Tully- 
Fisher relation we obtain in the Classic/SCDM models. The 
error bars indicate the 1-a variance over many Monte-Carlo 
realizations. In this plot, we have included only the model 
galaxies with more than 1O 7 M0 of cold gas, and which were 
classified as spirals a ccord ing to their bulge-to-disk ratio as 
described in Section 2.8.4. This is an attempt to select the 



model galaxies that most closely correspond to the galaxies 
in the observational Tully-Fisher samples we are consider- 
ing. The Munich and Santa Cruz packages with moderate 
feedback produce fairly good agreement with the slope and 
scatter of the observed TFR. Note that central galaxies tend 
to be brighter than the satellite galaxies. This is due to our 
assumption that all new cooling gas is accreted by the cen- 
tral galaxy, which may not be realistic. We intend to inves- 
tigate this using hydro simulations. The Durham package 
and the Santa Cruz package with high feedback both show 
curvature on the faint end due to the strong supernova feed- 
back. The curvature on the bright end of all of the Classic 
models occurs due to the static halo cooling model. 

Fig. |l| shows the TFR for the New/Santa Cruz (fidu- 
cial) models. The results are quite good for all of these mod- 
els. There is still a slight curvature on the bright end, but it 
is less pronounced, almost absent, in the low-f2 models. The 
models also show a bit of curvature at the very faint end, 
but this is beyond the level probed by the observations cur- 
rently under consideration. Comparison with samples that 
probe the TFR to fainter magnitudes is an important test 
of the supernovae feedback modelling. Note that the scatter 
also increases at fainter magnitudes, which is also observa- 
tionally testable. 

It should be kept in mind that this comparison rests on 
the assignment of model galaxy circular velocities as well as 
luminosities, and on the conversion from circular velocity to 
linewidth. In the current models we have assumed that all 
galaxies have perfectly flat rotation curves out to the virial 
radius of the halo; i.e., that the circular velocity measured 
by TF observations (typically at about two optical disk scale 
lengths) is the same as the virial velocity of the dark mat- 
ter halo. This assumption clearly must break down in halos 
with circular velocities larger than about 350-400 km s , 
as no known galaxies have rotation velocities this large. As 
we noted in Section 3, if the profiles nf Hark matter halos 



resemble the NFW profile, then V„ at a, few scale lengths will 



be larger than at r v ir for smaller (galaxy) mass halos, and 
smaller than at r v j r for larger (cluster) mass halos. In addi- 
tion, the dissip ative infall of baryons w i ll modify the inne r 



rotation curve (Blumenthal et al. 1986; Flores et al. 1993) 



Before attempting a rigorous evaluation of the Tully-Fisher 
relation predicted by the models, the halo and disk profiles 
should be modelled in more detail. We intend to address this 
problem in future work. 



10 



2 

CO 



d 



id 

a 

GO 




tCDM 




_L 



-16 



10 



03 

O 
O 



Q0 

O 



-16 -30 
M,-5 log h 



-33 



T 



ACDM.3 




W&f^:--^^: "■■■* 



*•>■:■ v. - 



-16 



-30 



-33 



M t — 5 log h 



Figure 16. Small dots show the cold gas masses of the model 
galaxies from mock catalogs extracted from the New/Santa Cruz 
(fiducial) tCDM (top) and ACDM.3 ("bottom ) models. The solid 



line is an approximate fit to Hi observations ( ie Blok, McGaugh 



& van der Hulst 1996). The dashed line is the same fit, with the 



gas masses multiplied by a factor of two to allow for a contribution 
from cold gas not in the form of H\. 



5.3 Cold Gas 

We now investigate the cold gas masses of galaxies in our 
models. This is an important counterpart to studying the 
luminosities of galaxies. Fig. |l^ shows the mass of cold gas 
in the model galaxies as a function of I magnitude for two 
examples of our fiducial models (rCDM and ACDM.3). The 
solid line shows an approximate fit to local Ht d ata (see 
Fig. 12 of|de Blok, McGaugh, fc van der Hulst 199rj). Recall 
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Figure 17. The mass function. The tw o smooth linns shnv 
the Schechter function fit to the results from Zwaan ct al. (1997 ) 



where the lower (bold) line is the actual fit and the upper line 
shows the effect of increasing all the gas masses by a factor of two 
to account for cold gas not in the form of Hj. The survey is sen- 
sitive in the range 10 7 to 10 10 A/q, and the dashed lines show the 
extrapolation of the Schechter function into the unprobed region. 
The thin bo ld lines show the u ncertainty in the faint end slope 
as given by [Zwaan et al. (1997 ) . The arrows show upper limits 
from a com plementary Arecibo survey, also mentioned in Zwaan 
et al. (l |997| ). The histograms show the results from the models. 
The top panel shows the Classic/SCDM models with different 
sf/fb packages, and the bottom panel shows the New/Santa Cruz 
(fiducial) models for different cosmologies. 
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Figure 18. The metallicity-luminosity relation in the Clas- 
sic/SCDM models (small dots), for galaxies within a "local group" 
(V c = 220 km s _1 ) sized halo. Bold lines show fits to the observed 
relation for bri ght spirals and for local dwa rf galaxies (from the 
compilation in Kobulnicky & Zaritsky 1998). 



that we set the free parameters to match the zero-point of 
this relation at Mr — 5 log h ~ —21.8, assuming that the mass 
of "cold gas" in our model reference galaxy is approximately 
a factor of two larger than the observed Hi mass to allow for 
molecular hydrogen (we neglect the additional contribution 
of helium and ionized hydrogen). This corresponds to the 
typical contribution of molecula r hydrogen in an Sb-Sc type 
galaxy (Young & Knezek 1989). The observations show a 
large scatter, comparable to the scatter in the models. The 
models results are consistent with the observed trend of gas 
mass with magnitude and the scatter in this relation. It 
should be noted that we have not made any morphological 
cuts on the model galaxies, whereas the observations are for 
late-type galaxies. The results look similar for all the models. 

We also investigate the Hi mass function, or the num- 
ber density of galaxies with a given Hi mass. Th is has been 
estimated by the survey of Rao fc Briggs (1993[). and more 



recent ly in the blind Hi survey described in 



Zwaan et al 



(1997). The latter should place strong upper limits on the 



number of low surface brightness galaxies (unless there is 
a very gas poor population) because it is not optically se- 
lected. In Fig. hTJ we show the Hi mas s function for all th e 



models, along with observations from Zwaan et al. (1997). 
All of the SCDM models show a considerable excess espe- 
cially on the small-mass end. The rCDM models show a 
somewhat smaller excess, and the other models are in good 
agreement with the observational limits across the range of 
Hi gas masses probed by the observations. 
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Figure 19. The mctallicity-luminosity relation for a volume lim- 
ited mock catalog extracted from the New/Santa Cruz (fiducial) 
ACDM.3 models. Large symbols show observations of bri ght spi 
rals, d warf irregulars a nd spheroidals, and HII galaxies (Kobul 
nicky fc|Zaritsky 1998|). 



5.4 Metallicity-Luminosity Relation 

Nearby galaxies are known to exhibit a trend between 
their B-band luminosities and their metal contents, in the 
sense that more luminous galaxies are more metal-rich. The 
slope of the observed relation derived for bright spirals 
{Mb — 5 log faso <: —18) is shallower than that for nearby 
dwarf galaxies ( gkillman. Kcnnicutt. fc Hodge 198S|; Richei 



& McC a.ll 19951; IZaritsky, Kcnnicutt, fc Huchra 1994 ^ oh 



urnicky ffcZaritsky 1998j ). 

Although we obtain a similar trend in the models (see 
Fig. [li]), the detailed behaviour of the observations is not 
well reproduced in any of the models. Recall that we have 
set our yield parameter y in order to obtain solar metallicity 
in our approximately "Milky Way" sized reference galaxy. 
The relation that we obtain depends on the treatment of 
metal and gas ejection by supernovae. In the Munich pack- 
age, none of the metal or gas is ejected from the halo, and 
this package produces a very shallow relation with a break 
at about Mb = —15. In the Durham package, all of the 
reheated gas and metals are ejected from the halo, result- 
ing in a very steep relation even for the bright galaxies. In 
the Santa Cruz package, the ejection of gas and metals is 
modelled using the disk-halo approach. This leads to a rela- 
tion which is consistent with the bright galaxies, but dwarf 
galaxies that are too metal-rich compared to the observa- 
tions. This is the case even in the high feedback package. 
The Durham package produces the best agreement with the 
observed relation. However, it also produced an unaccept- 
able degree of curvature on the faint end of the Tully-Fisher 
relation. 
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Figure 20. B-V colour histogram for galaxies in the New/Santa 
Cruz (fiducial) models, in the rCDM (top panel) and ACDM.3 
(bottom panel) cosmologies. The shaded histogram shows the 
observed colour distribution o f bright galaxies (Mr — 5 lo g h > 
-19.0) from the RC3 catalog (|De Vaucoulcurs ct al 199l[). The 



unshaded histograms show the model galaxies, selected to be 
brighter than M B - 5 log h = -19.0 (solid) or M B - 5 log h > 
—15.5 (dashed). Within each panel, the top half-panel shows the 
model results without the correction for dust reddening, and the 
bottom half-panel shows the results with the correction for dust. 



efficient than gas ejection. This has been proposed on the 
basis of high-resolution hydrodyna mic simulations of dwarf 
galaxies ( Mac Low fc Ferrara 1998 ) . Our results sup port the 



low the ejected gas may indicate that metal ejection is more 



strong ly differential ejection efficiency proposed by Martin 
(1998) on the basis of observations. One should use caution 
in interpreting the observations, however, because of pos- 
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sible systematic errors in the observational determination 
of metal abundances in different types of galaxies and us- 
ing diffe rent methods (Kobulnicky, Kennicutt, & Pizagno 
1999|)~~[nterpret ing the observations is also complicated by 
the presence of metallicity gradients in large galaxies, and 
t he strong cor relation of metallicity with surface brightness 
((Harriett 1997). We postpone a more careful investigation of 



these issues to future papers. 



5.5 Colours 



A familiar property of observed galaxies is the "colour- 
magnitude" relation: bright galaxies are observed to be red- 
der than fainter ones. It is an often repeated statement that 
hierarchical models of galaxy formation generically predict 
that more massive objects form "later" than smaller mass 
objects. This statement is often misinterpreted to imply that 
larger mass objects should be "younger" and therefore bluer 
than smaller ones. If the formation time of an object is de- 
fined as the time when a given fraction of its mass has been 
assembled into a single progenitor, then it is true that larger 
mass objects have later formation times than smaller mass 
objects. However, if we define "age" as the time spent in 
the sort of environment where we expect that star forma- 
tion is able to occur (i.e., within a collapsed halo), then the 
mean age of the material in large mass halos is older than 
in smaller mass halos. This is because large mass halos are 
associated with higher peaks in the density field, which col- 
lapse earlier. 

In the SAMs, we find that when the effects of dust and 
metallicity are neglected, we obtain flat colour-magnitude 
relations or (depending on the model and the colour bands 
in question) bright galaxies that are only slightly redder (see 
Fig. top panel) than the faint population. The bright 
galaxies are also a bit too blue overall compared with obser- 
vations (the obser ved colour distribution of brig ht galaxies in 
the RC3 catalog (De Vaucouleurs et al. 1991) is shown for 
comparison). Inclusion of dust extinction shifts the colour 
distribution towards the red, and shifts the bright galaxies 
more than the faint ones because of the differential nature 
of our dust recipe. This brings the optical colours int o fairly 



good agreement with the observations. In addition, Kauff- 



inn & Chariot (1996) have shown that the inclusion of 
metallicity effects on the model spectra (here we have used 
only the solar metallicity stellar population models) can also 
produce the observed colour-magnitude slope. Presumably 
the observed trend is a combination of these two effects. 



5.6 Sizes 

We estimate the exponential scale radii of ou r model disks 
using the approach described in Section 2.4. We show the 
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relatior ship between scale radius and circular velocity in 



Fig. 21 The upper left-most panel shows the observations 
of pourteau (1996[ ) for late-type galaxies, where V c is the 
disk rotation velocity at 2.2 scale-lengths, and r„ is the ex- 
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Figure 21. Exponential scale radius vs. circular velocity. The 
large filled squa res (top left panel) show the observations of 
Courteau (199tj ). The small dots show the results of the fidu- 
cial Santa Cruz models in various cosmologies as indicated. Bold 
solid lines show the relation we would obtain if all gas fell in from 
the virial radius of the halo, for Ah = 0.05 (50 percent point). 
Bold dashed lines show the same relation for Ah = 0.025 (10 per- 
cent) and Ah = 0.1 (90 percent). These vary slightly depending 
on cosmology. The rCDM (lower) and ACDM.3 (upper) relations 
are shown on the upper left panel with the observations as a ref- 
erence point (other models are intermediate). The arrows on the 
rCDM and ACDM.3 panels show the corrections to this relation 
predicted from the fitting formulae of Mo, Mao, & White (1998) 
(see text). 



may be smaller yet. The bold solid lines indicate the relation 
r s — l/\/2A£fr v i r , which we would obtain if the gas fell in 
from the virial radius of the halo. As we have discussed in 
Section 



2.4 



we have used a constant value of Ah = 0.05 for 
all halos, so the only source of scatter in the relation that we 
obtain is from the different cooling radii of the halos, which 
occur as a result of the scatter in halo merging histories (see 
Figure [J). The SAM dots always lie at smaller r s than the 
bold lines because the cooling radius is always smaller than 
the virial radius. The bold dashed lines show the same rela- 
tion for Ah = 0.1 and Ah = 0.025, which are the 10 and 90 



percent points of the distribution of Ah used by Mo, Mao 



& White (1998) 



As we discussed in Section p.4 
"back reaction 



we have not included the 



of baryons on the dark matter during their 
collapse, which will tend to lead to smaller scale radii (|Blu 



al scale length in the r band. In the models, the scale 


menthal et al. 1986 




Flores et al. 1993 




Mo, Mao, & White 


hat we estimate represents the total baryonic mass 


1998 ) and will also modify the rotation 


curve. The use of a 



(stars and cold gas) in the disk. The scale radius of the stel- 
lar mass may be smaller if star formation is more efficient in 
the inner parts of the disk, and the scale of the optical light 



more realistic halo profile (e.g. NFW) will also change these 
results. The correction due to these effects, as pr edicted by 



the fitting formulae of Mo, Mao, fc White (1998), is shown 
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by the arrows in the rCDM and ACDM.3 panels (the cor- 
rection is relatively insensitive to cosmology, and is similar 
for the other models). The direction of the correction is to 
produce galaxies with larger V ma x and smaller scale radii. 
The details of these corrections depend on the assumed halo 
profile and disk baryon fraction as well as other parame- 
ters. We intend to incorporate improved modelling of disk 
sizes and rotation curves into our models and present more 
detailed predictions of these quantities in the near future. 

In the meantime, several things are worth noting. An 
obvious difference in the current model predictions is the 
break in the r a -V c relation at about 200 km s _1 : the trend 
reverses and the scale radii start to decrease at larger V c . 
This is caused by the decreasing cooling radius in halos with 
higher virial temperature (see Fig. ^). No such break is ev- 
ident in the observations, although the sample contains a 
relatively small number of galaxies with V c > 200 km s _1 . 
The rightward shift indicated by the arrows (due to "peak- 
ing up" of the rotation curve caused by the effects mentioned 
above) may solve this problem, but will also make the disks 
too small at a given V c compared to the Courteau data. 
Most of the effects we have mentioned indicate that we may 
already be systematically over-estimating the disk sizes. We 
therefore may be facing a puzzle similar to the "angular 
momentum" problem found in N-body simulations with hy- 
drodynamics (i.e., disks are too small and concentrate d at a 



given c ircular velocity compared to observations, cf. stein- 
metz 



Navarro 199S) 



6 SUMMARY AND DISCUSSION 

We have presented new semi-analytic models of galaxy for- 
mation and shown that these models can reproduce many 
key observational properties of galaxies in the lo cal Universe. 



Our approach is similar to that i ntroduced by Kauffmann 
White, fc Guiderdoni (199^ ) and |Cole et al. (1994| ), but 



have introduced several modified or new ingredients, includ- 
ing: 

• Somerville-Kolatt met hod for "planting" merger trees 
(somerville & Kolatt 1998; shown to give good agreement 
with merger trees ex tracted from N-body simulations in 
Somerville et al. 1998|). 



• Improved Sheth-Tormen model ( 3heth & Tormen 199£ ) 
for the mass function of dark matter halos (an improved ver- 
sion of the Press-Schechter model which gives much better 
agreement with the reusults of N-body simulations). 

• "Dynamic halo" model for gas cooling (includes the ef- 
fects of halo merger events on the density and temperature 
of the hot halo gas). 

• "Disk-halo" model for supernovae feedback (models the 
ejection of cold gas from the disk and global (dark matter) 
potential seperately) 



• Dust extinctio n based on the empirical recipe of Wang 
fc Heckman (1996| ) 

• Galaxy mergers due to satell ite collisions, using t he 
simulation-based approximation of Vlakino & Hut (1997) 



• More detail model ling of starbursts based o n hydro- 
dynamical simu lations (Mihos & Hernquist 1994; Mihos & 
Hernquist 1996b 



We have investigated several different "packages" of recipes 
for star formation and supernovae feedback in order to gain 
better understanding of the importance of the way in which 
these processes are parameterized. We have also illustrated 
the results of varying the most important free parameters in 
our models. 

We have addressed the long-standing problem of the 
physical explanation of the observed luminosity function 
within CDM-type hierarchial models of structure forma- 
tion. Early in the history of CDM, it was noted that the 
mass function of dark matter halos, whether predicted by 
analytic models like Press-Schechter or derived from N- 
body simulations, has a steep power-law slope (a ~ —2) 
for masses < 10 14 /i _1 Mq, and an exponential cut-off at 
~ lO 14 /!" 1 Mq, much larger than the expected mass of the 
halos surrounding L„ galaxies, as estimated from their inter- 
nal velocity dispersions. It was proposed that feedback due 
to supernovae could suppress star formation in small mass 
halos, leading to a flatter faint end slope, and inefficient gas 
cooling in large mass halos could cause the "knee" at £*. 
However, the first generation of SAMs, which attempted to 
actually model these processes in some detail, encountered 
some difficulties. The Munich models produced the correct 
(B-band) TFR slope, but the faint-end slope of the luminos- 
ity function was still too steep. In addition, unless an ad-hoc 
cutoff was applied, in which gas cooling was turned off by 
hand in halos larger than 500 km s , these models did not 
produce a "knee" in the luminosity function and showed an 
excess of very bright galaxies. These models were normal- 
ized to the observed luminosity of the Milky Way Galaxy, 
and were claimed to reproduce the observed zeropoint of the 
local B-band TFR. The Durham models produced a lumi- 
nosity function with a "knee" , and the free parameters were 
adjusted in order to match its location with that of the ob- 
served B-band luminosity function. Their luminosity func- 
tions showed a flatter faint-end slope, in better agreement 
with observations, but produced a TFR with a zeropoint 
offset of about 2 magnitudes and a serious deviation from 
the observed power-law behavior on the small V c end. 

We have clarified the reasons for some of these differ- 
ences. First, we have explained how it is that, although the 
Munich group claimed to reproduce the observed zeropoint 
of the TFR, in fact their model galaxies were ~ 2 magni- 
tudes too faint at a given circular velocity compared to the 
TFR derived from recent large I-band samples. Using our 
models with sf/fb recipes chosen to be similar to those of 
the Munich and Durham group, we showed that when the 
models are normalized in the same way, the luminosity func- 
tion and TFR are nearly identical for bright/large V c galax- 
ies (M B - 5 log h <; -20, V c £ 220 km s" 1 ). The results 
differ substantially only for faint/small V c galaxies, and this 
difference can be traced mainly to the stronger supernovae 
feedback recipe assumed by the Durham group. We demon- 
strate this by showing that as we turn up the parameter 
that represents the fraction of supernovae energy deposited 
in the cold gas [e% N ), our results move continuously from a 
situation resembling the Munich models (steep LF, power- 
law TFR) to one resembling the Durham models (flat LF, 
curved TFR). This works because in our "disk-halo" feed- 
back model, the parameter e° SN affects not only the total 
amount of gas that is reheated, but also the fraction that 
is ejected from the halo. We issue a warning, however, that 
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although the supernovae feedback efficiency is the dominant 
factor determining these results at z = 0, the redshift evo- 
lution is also very sensitive to the assumed star formation 
recip e. This will be illustrated in detail in a companion pa- 
per (Somerville, Primack, & Faber 1999). 

We have shown that the observed B and K band lumi- 
nosity functions and the Tully-Fisher relation can be repro- 
duced simultaneously in our models. This improvement is 
not due to any one effect but is the result of many combined 
factors, as summarized below. We first normalize our models 
to fix the zeropoint of the TFR for a typical L* galaxy. This 
requires us to use a stellar mass-to-light ratio approximately 
a factor 2 to 3 times higher than the published Munich and 
Durham models, but corresponds to taking the predictions 
of the Bruzual- Chariot stellar population models at close to 
face value and is in better agreement with observational es- 
timates. If we then use cooling, feedback, and star formation 
recipes similar to the published Munich models, we obtain a 
B-band luminosity function with several problems: the over- 
all normalization is too high, the faint end slope is too steep, 
and there is a "tail" of bright galaxies (see Fig. model 1). 
Fig. ^ shows the effect of introducing various changes, one 
by one, which eventually lead to our fiducial model choices. 
These are summarized below. We remark whether the effect 
is important on the bright (Mb — 51og/i ~ —21.5) or faint 
(Mb — 5 log h ~ —16.5) end of the luminosity function, and 
by what factor the B-band luminosity function changes at 
this magnitude. These factors should be considered approx- 
imate only, and may be read from the top panel of Fig. ^2| 

• "disk-halo" feedback model (faint end, factor of 2.5) 

• Sheth-Tormen mass function (overall, factor of 1.5) 

• "dynamic halo" cooling model (bright end, factor of 3) 

• low Q ~ 0.5 cosmology (bright end, factor of 1.6; faint 
end, factor of 3) 

• dust extinction (bright end, factor of 16) 

Note that the observations have also changed — the solid 
bold curve in Fig. [22] shows t he luminosity functi on derived 
from the recent 2dF survey jFolkcs et al. 1999| ), which is 
i n excellent agreem ent with the LF from the deep ESO slice 
( Zucca et al. 1997| ). Both have considerably steeper faint-end 



' B— band 



slopes than the LF derived from the APM survey Loveday 
et al. ( 992), which was the standard at the time of much 



of the earlier modelling. The more recent observations are 
easier to reconcile with the models. A similar accounting 
may be done for the K-band luminosity function (Fig. |22l 
bottom panel). 

We therefore conclude that the very strong feedback and 
suppression of star formation in small V c galaxies assumed 
in the Durham models is not necessary in order to repro- 
duce the observed luminosity function, and is disfavored as 
it produces curvature on the small V c end of the TFR. 

In our n = 1 models (SCDM and tCDM), we find that 
in order to produce galaxies with large enough luminosi- 
ties and gas masses, we must assume values of the baryon 
fraction (/baryon ~ 0.1 — 0.12) which are rather high com- 
pared w ith estimates from observations of high-redshift deu- 
terium ( Tytler et al. 1999 ) , though consistent with estimated 



baryon fractions in groups and clusters. In our best fio = 1 
models, we find good agreement with the general shape of 
the LF, but the model LF is too high by an overall factor 
of ~ 3 in B and 2.5 in K. The mass function of cold gas is 
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Figure 22. The effect of various model variations on the model B 
and K-band luminosity functions, introduced one by one. The new 
feature of each model is in italics. Model 1: Munich-style feed- 
back, original Press-Schechter mass function, static halo cooling, 
SCDM, no dust correction. Model 2: disk-halo feedback, orig- 
inal Press-Schechter mass function, static halo cooling, SCDM, 
no dust correction. Model 3: disk-halo feedback, Sheth-Tormen 
mass function, static halo cooling, SCDM, no dust correction. 
Model 4: disk-halo feedback, Sheth-Tormen mass function, dy- 
namic halo cooling, SCDM, no dust correction. Model 5: disk- 
halo feedback, Sheth-Tormen mass function, dynamic halo cool- 
ing, ACDM.5, no dust correction. Fiducial model: disk-halo 
feedback, Sheth-Tormen mass function, dynamic halo cooling, 
ACDM.5, with dust correction. In the top panel, the bold lines 
show the observed B-band luminosity function from the 2dF sur- 
vey (higher line) and the APM survey (lower line). In the bottom 
panel, the bold lines show the obser ved K-band luminosi ty func- 
tion from Szokoly et al. (199gj ) and Gardner et al. (19971). 
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also a factor of ~ 5 higher than estimates from blind Hi sur- 
veys. In order to reconcile these models with observations, 
we would have to believe that there is a substantial popu- 
lation of galaxies, including some with large total masses, 
that are undetected in optical emission or radio emission 
from cold Hi gas. This seems unlikely, though not impossi- 
ble. Thus, although we cannot say that Qq = 1 is ruled out, 
our results are certainly more easily compatible with mod- 
els in which Qq ~ 0.3 — 0.5, with or without a cosmological 
constant. 

The same fiducial models produce good agreement with 
observations of the mass function of cold Hi gas and the 
magnitude-Hi-mass relation. When we normalize our models 
to produce a "Milky- Way" galaxy with solar metallicity, the 
metallicities of dwarf galaxies in our models are somewhat 
higher than the average metallicity of nearby dwarf galax- 
ies, i.e. the slope of the metallicity- luminosity relation is too 
shallow. This may be evidence that metals are ejected by su- 
pernovae more efficiently than the cold gas, or an indication 
that our "constant yield" approach to modelling chemical 
evolution is too simplistic. Alternatively, it may be due to 
systematic uncertainties in deriving observational es timates 



of met al abundances in different types of galaxies ( Kobul- 
nicky, f ^ennicutt, fc Pizagno 199£ ). Our fiducial models pro- 



duce good qualitative agreement with the optical colors of 
bright galaxies, and reproduce the observed color-magnitude 
trend, when dust extinction is included. Although the rela- 
tionship between the exponential scale radius and circular 
velocity of disks that we estimate is in reasonably good qual- 
itative agreement with observations, we conclude that more 
detailed modelling is necessary. 

A great strength of the SAM technique is that one can 
make self-consistent predictions pertaining to a wide variety 
of observations. In companion papers, we investigate our 
predictions for the properties of high-redshift galaxies and 
t he history of stars, cold gas, and m etals at high redshift 
( Somerville, Primack, fc Fabcr 1999|) , and extend our pre- 
dictions for local galaxies to shorter (far UV) and longer (far 



IR to s ub-mm) wavelengths ( Somerville et al. 1999| ; 



Somerv ille, fc Primack 199< ) 



Bullock 
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Figure Al. The relationship between halo mass and virial ve- 
locity from the spherical tophat model, at z = (bottom set of 
lines), 2 = 1, and 2 = 3 (top), for the cosmologies discussed in the 
text. The relation depends (weakly) on cosmology and (strongly) 
on redshift. 



APPENDIX A: SPHERICAL COLLAPSE IN A 
GENERAL COSMOLOGY 

We need to be able to relate the mass, radius, and veloc- 
ity dispersion of dark matter halos for any given redshift. 
This is made possible using the spherical collapse model, 
one of the apparently gross oversimplifications that seems 
to work surprisingly well. We imagine a spherical patch of 
the universe with a uniform overdensity Si within a radius 
n at a very early time tj (often called a "top-hat" perturba- 
tion). We assume that the collapsing shells of matter do not 
cross. If we consider a particle at radius r, Birkhoff's the- 
orem ( |Birkhoff 1923TJ tells us that we can ignore the mass 
outside this radius in computing the motion of the particle. 
The equation of motion for our particle (in physical, rather 
than comoving, coordinates) is then 



GM 



A 

+ r 



(Al) 



where M = (47r/3)rf /0& (ti) (1 + Si) and pbiti) is the back- 
ground density of the universe at t,. Integrating this equa- 
tion gives 



r = H 



^(l + «y% + fi A r 2 



K 



(A2) 



where K is a constant of integration. We may fix this by 
noting that if we have picked ti early enough that Q ~ 1 at 
that time, linear theory tells us that the initial velocity is 



r(U) = H, 



—3- H 5- + "A- 



(A3) 



(Peebles 1984). At the point of maximum expansion, or 
"turnaround", r — 0. If we set equation (A3) to zero, we 
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obtain a cubic equation for rt a , the radius of the perturba- 
tion at turnaround, which must be solved numerically for 
the general cosmology given here, but for specia l cases it 
can be solved analytically (cf. Padmanabhan 1993| ), From a 
symmetry argument, we note that the time when the per- 
turbation collapses to a point, i co ii, is always twice tta (the 
time at maximum expansion). We can now write an implicit 
equation for the mass of a perturbation that is collapsing at 
icon : 



icon — 2 / — . 
Jo 



r 



(A4) 



We know the mass and the radius at turnaround, so we can 
calculate the density of the perturbation at turnaround, pta- 
Of course the perturbation will not really collapse to 
a point. Before that happens, shell crossing will occur, and 
it will virialize. We can find the radius after viralization in 
terms of the turnaround radius using the virial th eorem. The 
total energy at turnaround is (Lahav et al. 1991) 

r, „ rr 3 GM 2 
E = U G ,ta + U A ,ta = ~- 

5 r ta 

where the second term is due to the cosmological constant. 
Now using the virial theorem for the final state: 



- To AM ^ 



(A5) 



T 



U G j + U A j ■ 



(A6) 



From conservation of energy we then have \Ugj + 2J7a,/ = 
Ua,ta + UA,ta- This leads to a cubic equation for the ratio of 
the virial radius r v i r to the turnaround radius r ta - We now 
know r v ir and can write down the virial density 



A c (z) 



gviffl(z) 

n pg(i + zY ' 



(A7) 



We now have a relationship between the mass, virial 
radius, and collapse redshift z. If we assume a radial profile 
for the virialized halo, we can use the virial theorem again 
to relate these quantities to the velocity dispersion. If we 
assume that the halo is a singular isothermal sphere, p oc 
r~ 2 , truncated at the virial radius, then we have 



J 2 

— a = 
2 



GM Arl 



2r vir 18 

or, in terms of the circular velocity V c , assuming V 2 

GM n A 



Vc 



. u-2 2 

-H r vh 



(A8) 



2a 2 



(A9) 



We can now translate between mass and velocity dispersion 
at any given redshift. Note that in universes with a non-zero 
cosmological constant, halos of a given circular velocity are 
less massive because of the A contribution to the energy. 



In pra ctice, we use the fitting formula of Bryan & Nor- 
man (1 )97 ) for the virial density: 



A c = 18tt + 82x - 39z 



for a flat universe and 
A 



18tt 2 + 60:r - 32a; 2 



(A10) 



(All) 



for an open universe, where x = S}(z) — 1. This formula is 
accurate to 1% in the range 0.1 < fi < 1, which is more 
than adequate for our purposes. We now can write down 
the general expression for r v i r in closed form. 



M 



1 1/3 



47T A c (z)n pc,0 



(A12) 



In conjunction with equation A9, this allows us to calculate 
the circular velocity and viral radius for a halo with a given 
mass at any redshift z. These expressions are valid for open 
cosmologies with A = and flat cosmologies with non-zero 
A. 
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